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PREFACE 


For  over  8 years , the  author  has  conducted  a continuing 
research  project  that  has  produced  new  computer  progru*s , some 
of  which  have  wide  applicability  to  allocation  of  resources 
(such  as  weapons)  to  objectives  (such  as  targets).  Some  of 
the  programs  are  general  in  nature  and  have  been  used  by 
other  ANSER  analysts  as  well  as  the  author  to  solve  a variety 
of  weapon  systems  analysis  problems  in  our  work  for  the 
Deputy  chief  of  Staff/Research  anu  Development,  Headquarters, 
United  States  Air  Force.  Several  programs  have  also  been 
given  to  various  government  agencies  or  defense  contractors. 
Documentation  of  these  has  been  minimal. 

This  report  describes  companion  programs,  PATH87  and 
PATH87A,  that  are  very  general  in  nature  and  which  are 
milestones  in  a long  evolutionary  development.  This  report 
was  originally  drafted  under  the  same  title  in  March  1971, 
but  the  final  version  was  not  completed  because  of  the 
press  of  other  business.  However,  it  was  circulated  in 
draft  form  to  a number  of  defense  contractors,  some  of  whom 
were  kind  enough  to  list  the  draft  as  a reference  in  their 
own  publications.  The  current  edition  represents  a minor 
revision  and  update  of  the  1971  draft. 

It  has  been  demonstrated  that  the  PATH  method  offers 
unique  advantages  of  speed  and  flexibility  in  solving 
problems  facing  defense  analysts  and  it  is  hoped  that 
publication  of  this  report  through  the  National  Technical 
Information  Service  of  the  Defense  Documentation  Center 
will  make  this  method  more  widely  available.  Also,  the 
method  has  features  which  can  be  applieu  to  many  problems 
of  resource  allocation  facing  nondefense  planners. 
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MULTIDIMENSIONAL  PREFERENTIAL  STRATEGIES 


I.  INTRODUCTION 

This  report  presents  a method  for  solving  a large  class 
of  resource  allocation  problems  that  are  multidimensional  in 
the  sense  ot  having  many  types  of  resources  and  many  types  of 
objects  to  which  resources  are  allocated.  The  class  of 
problems  is  restricted  to  those  in  which  discrete  quantities 
of  resources  are  allocated  to  any  object,  for  example,  inte- 
gral numbers  of  weapons  to  a target.  The  solution  method 
involves  parametric  linear  programming. 

Historically,  the  method  was  developed  in  order  to 
generalize  earlier  ANSER  ^ork  (References  1 and  2) . The 
essential  features  ot  Sections  II  and  III  were  the  subject 
of  an  oral  seminar  at  tho  ORSA  meeting,  1-3  May  1968,  San 
Francisco,  California.  (Paper  TP  1.12,  A Generalized  Weapon 
Allocation  Game.) 

A.  The  Problem 

The  general  problem  of  the  class  with  which  we  are  deal- 
ing is  formulated  as  a two-sided  mathematical  game  in  Section 
II.  The  formulation  includes  one-sided  optimizations  as 
special  cases  where  one  of  the  players  of  the  game  has  zero 
resources . 

An  important  feature  of  the  formulation  is  the  definition 
of  an  allocation,  or  strategy,  in  such  a way  that  the  solution 
variables  appear  linearly  with  constant  coefficients  in  the 
objective  function  <=nd  the  constraints.  Any  nonlinear  features 
of  the  problem  arc  incorporated  in  the  computation  of  those 
coefficients . 

Physical  feasibility  of  solutions  is  ignored,  and  the 
solution  strategies  may  contajn  fractional  values  that  are 
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apparently  not  feasible.  However,  these  fractional  strategies 
can  be  realized  by  one  or  more  precisely  equivalent,  feasible, 
mixed  strategies  in  nearly  every  two-sided  game  of  practical 
interest.  In  most  one-sided  optimizations,  the  physical 
infeasibility,  if  any,  is  a local  one  of  no  great  significance 

Solution  strategies  are  obtained  by  solving  simultaneous 
linear  equations,  with  the  matrix  of  coefficients  for  one 
player  being  the  transpose  of  that  for  the  other  player. 

This  common  matrix  constitutes  a basis  tor  the  solution, 
called  here  a Q-b.  sis. 

B.  Path  Method 

The  use  of  parametric  linear  programming  to  determine  the 
Q-basis  is  discussed  in  Section  III. 

Beginning  with  some  set  of  resources,  usually  zero,  for 
which  the  Q-basis  is  known,  resources  are  then  varied  con- 
tinuously in  any  arbitrarily  prescribed  manner;  that  is,  along 
any  prescribed  path  in  the  space  of  resources.  The  path 
method  makes  appropriate  changes  in  the  Q-basis  at  critical 
points  of  the  path. 

Solutions  for  the  two  players  are  alternated  along  the 
path,  with  a resource  parameter  being  used  for  tha  player 
whose  resources  are  changing  and  a marginal-value  parameter 
for  the  other  player. 

C.  Computer  Programs 

Two  computer  programs  are  included  in  the  report.  The 
PATH87  program,  designed  to  solve  multidimensional  two-sided 
weapon-allocation  games  for  point  targets,  is  listed  and 
discussed  in  Section  IV.  This  program  is  the  1971  version 
in  an  evolutionary  sequence  pointed  toward  greater  capacity, 
speed,  flexibility,  and  reliability. 
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The  other  program  is  an  abbreviated  version  designed  to 
optimize  one-sided  weapon  allocations  for  point  targets. 
Called  PATH87A,  it  is  listed  and  discussed  in  Section  V. 

Readers  who  are  not  interested  in  the  details  of  computer 
programming  should  skip  Sections  IV  and  V. 

D.  Preferential  Strategies 

Illustrative  examples  of  actual  runs  of  both  programs 
are  given  in  Section  VI.  One  purpose  of  the  section  is  to 
show  someone  who  may  have  skipped  the  earlier  sections  how 
to  use  the  programs,  how  to  input  data,  and  how  to  read  the 
output.  A second  purpose  is  to  point  out  and  explain  salient 
characteristics  of  typical  solutions. 

E.  Other  Applications 

The  path  method  applies  to  a variety  of  resource  alloca- 
tion problems  that  can  be  solved  by  suitable  modifications 
of  the  two  programs.  Section  VII  lists  some  features  of 
problems  that  have  been  solved  and  discusses  the  modifica- 
tions in  general  terms. 
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II.  THE  PROBLEM 


A multidimensional  resource  allocation  problem  of  the 
class  with  which  we  are  concerned  can  be  formulated  as  a 
matrix  game.  In  this  section,  we  will  give  a general  formu- 
lation of  the  game,  including  general  solution  conditions. 

We  will  then  develop  more  detailed  conditions  in  the  form  of 
a Q-basis,  consisting  of  a certain  coefficient  matrix  and 
some  associated  vectors.  Finally,  we  will  show  that  a 
solution  can  be  computed  directly  from  the  Q-basis  if  cer- 
tain conditions  are  satisfied. 

A.  Problem  Formulation 

1 . Players,  Resources,  and  Objects 

In  general  terms,  there  are  tv/o  opposing  players,  each 
having  resources  of  different  types  to  be  allocated  to 
specific  objects  (or  activities)  of  which  there  are  differ- 
ent types.  One  player  is  called  the  minimizing  player  because 
he  seeks  an  allocation  that  will  minimize  some  common  measure 
of  value.  His  opponent  is  called  the  maximizing  player. 

Each  player  must  make  his  own  allocation  in  ignorance  of  his 
opponent's  specific  allocation,  but  he  does  know  all  the 
resources  available  to  the  other. 

For  illustration,  an  attacker  would  try  to  allocate  a 
variety  of  weapons  to  a variety  of  targets  in  such  a way  as 
to  minimize  the  expected  value  of  the  surviving  targets, 
whereas  a defender  would  try  to  allocate  defensive  weapons 
so  as  to  maximize  expected  surviving  value. 

Let  us  adopt  the  following  notations 

M = number  of  types  of  minimizing  resources 

N « number  of  types  of  maximizing  resources 

Preceding  page  blank 
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G - number  of  types  of  objects. 

Let  us  also  represent  the  number  of  units  of  each  type  by: 

A„  for  m s 1 , • • • , M 
m 

Dn  for  n = 1,  • • • , N 

Tg  for  g ■ 1 , • • * , G 

Most  of  our  subsequent  analysis  is  concerned  with  the 
general  case  where  both  players  have  some  resources.  How- 
ever, one  should  observe  that  the  analysis  is  also  valid  if 
one  of  the  players  has  no  resources,  in  which  case  the  problem 
becomes  a one-sided  optimization.  As  a convention,  we  will 
suppose  this  case  to  be  represented  by  N = 0,  so  that  we 
will  always  have  M > 1 and  G > 1. 

2 . Elementary  Strategies 

An  elementary  strategy  is  defined  as  any  allocation  by 
one  player  to  a single  object  of  a particular  type.  It 
consists  of  some  number  of  units  of  each  of  that  player's 
resource  types.  We  will  denote  an  elementary  strategy  for 
the  minimizing  player  by  the  vector. 


5 (ail' 


*L' 
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and  one  for  the  maximizing  player  by  the  vector, 


<5?  = (<$?.  , 
3 1] 


6g. 

»3 


6g.) 


The  superscript  g identifies  the  object  type  to  which  the 
elementary  strategy  applies.  The  subscript  i (or  j)  identi- 
fies the  particular  elementary  strategy  among  all  those 
applying  to  object  type  g,  under  the  assumption  that  they 
are  arranged  in  some  numerical  order,  the  actual  arrangement 
being  immaterial.  The  subscript  m (or  n)  distinguishes  the 
resource  type,  so  that  the  component  is  the  number  of  type 
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resources  in  the  i elementary  minimizing  strategy  for  object 
type  g.  We  will  call  an  elementary  strategy  with  all  zero 
components  a null  strategy. 

3.  General  Strategies 

We  will  now  define  a general  strategy  (or  simply  a 
strategy)  for  either  player  as  a specific  allocation  of  all 
his  resources  among  all  the  objects.  A strategy  can  be 
described  in  terms  of  elementary  strategies  by  specifying 
the  fraction  of  the  number  of  objects  of  each  type  on  which 
each  elementary  strategy  is  used. 

For  the  minimizing  player,  we  will  denote  by  x?  the 
fraction  of  the  number  of  objects  of  type  g on  which  the 
elementary  strategy  u?  is  used.  The  components,  x?,  must 
satisfy  three  conditions:  every  component  must  be  non- 

negative; every  object  must  have  some  elementary  strategy, 
if  only  the  null  strategy;  and  all  resources  must  be  used. 

These  conditions  can  be  represented  by: 

x?  1 0 for  every  g,i 

Z x?  = 1 for  every  g (1) 

i 1 

Z Z T a?  x?  = A for  every  m 
. g im  i m 

g i ^ 

If  the  components  are  arranged  sequentially,  they  form  a 
s tr a tegy  vector: 

X = (x1,  x1,  *•♦,  •••,  x?,  •••)  , 

12  A ^ 

which  defines  an  allocation  of  minimizing  resources  to  objects. 

In  similar  fashion,  a strategy  vector  for  the  maximizing 
player  is  given  by: 

Q 0 

Y = (yj,  y\,  *•*,  y?»  Yy  *••)  * 
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with  conditions/ 


y?  i 0 for  every  g,j 

Z yg  - 1 for  every  g (2) 

j 3 

l I T 6g . yg  = D for  every  n 
g j 9 nj  *3  n 

At  this  point,  we  may  observe  that,  in  a one-sided  optimization 
problem  (characterized  by  N = 0) , the  only  allowable  elementary 
maximizing  strategy  on  each  object  is  the  null  strategy, 

6g  - (0) . With  this  restriction,  yg  = 1,  and  the  only  maxi- 
mizing strategy  vector  becomes 

G times 

Y » (1,  ••  •,  1)  . 


4 . Value  Functions 

The  statement  of  a problem  must  provide  some  way  for 
evaluating  the  worth  of  a strategy  and  making  a choice  of 
the  "best"  strategy.  The  basis  of  this  evaluation  is  the 
value  function,  one  of  which  must  be  specified  for  each 
object  type.  Each  value  function  must  define  a value 
associated  with  every  allowable  combination  of  an  elementary 
minimizing  strategy  and  an  elementary  maximizing  strategy  on 
tho  same  single  object.  Thus,  if  elementary  strategies 
u?  and  6?  occur  on  the  same  object  of  type  g,  the  value  may 
be  defined  as  v? ^ . The  value  function  for  object  type  g 
then  consists  of  the  matrix, 

vg  = (v? . ) 

There  must  be  a different  value  function  for  each  object 
type;  otherwise  object  types  could  be  combined  to  reduce  the 
dimensions  of  the  problem.  The  several  value  functions  must 
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be  expressed  in  terms  of  a common  unit/  or  measure  of  value, 
but  they  need  not  have  similar  formulations.  Indeed,  an 
explicit  formulation  is  not  required  at  all,  and  a table 
of  arbitrary  values  will  serve  the  purpose. 

There  is  great  latitude  for  specifying  value  functions. 

For  example,  in  a weapon  allocation  problem,  an  initial  value 
might  be  specified  for  a single  target  of  each  type  and  a 
formula  given  for  computing  expected  surviving  value  for  any 
pair  of  elementary  strategies.  The  value  function  in  a 
transportation  problem  might  be  the  cost  matrix  of  mileages 
between  sources  and  destinations. 

Figure  1 illustrates  the  relationship  between  elements  of 
a value  function  and  associated  strategies  for  a single  object 
type.  Figure  2 illustrates  a useful  concept  of  the  combined 
value  functions  and  strategies,  where  the  blocks  for  each 
group  symbolize  sets  of  elements  like  those  in  Figure  1. 


5.  Objective  Function 

When  the  opposing  players  use  specific  general  strategies, 
X and  Y , each  in  ignorance  of  the  other's  strategy,  there  is 
an  aggregate  expected  value  which  is  a function  of  the 
strategies  and  is  called  the  objective  function. 


In  developing  the  objective  function,  let  us  consider 
first  a single  object  of  type  g,  and  suppose  a?  is  the 
elementary  minimizing  strategy  on  this  object.  The  probabil- 
ity that  the  elementary  maximizing  strategy,  6?,  is  used  on 
this  object  is  y? , i.e.,  the  fraction  of  objects  of  this  type 
on  which  the  strategy  <5?  is  used.  The  expected  value 
associated  with  this  single  object  is  then  given  by 


E 
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FIGURE  2 

COMBINED  VALUE  FUNCTIONS 
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The  number  of  objects  on  which  the  strategy  a?  is  used  is 
Tg  xg.  Hence,  the  sum  of  the  expected  values  on  all  objects 
of  typ^  g is  given  by 


IT  xg  L v? , y? 
i g 1 j x3  3 

The  sum  of  these  values  over  all  object  types  is  the  objective 
function,  which  can  be  written 

V(X,Y)  ■ 1 I T.  x?  (T  v?.)  yg  . 

g i j 1 g *3  3 

This  is  the  function  that  one  player  seeks  to  minimize  by 
choosing  a strategy  X subject  to  the  constraints  (1).  The 
other  player  seeks  to  maximize  i t by  choosing  a strategy  Y 
subject  to  the  constraints  (2) . The  problem  of  finding 
best  strategies  for  each  player  is  a mathematical  game. 

In  the  foregoing  formulation,  the  objective  function  is 
linear  in  the  components  of  either  the  X or  Y strategy,  which 
permits  linear  programming  methods  to  be  used  for  the  solution. 
A similar  formulation  for  one-sided  optimizations  has  been 
used  by  Dianich  and  Hennig  (Reference  3) . 

6.  Solution 

For  a solution,  we  adopt  the  standard  criterion  that  a 
best  strategy  for  the  minimizing  player  is  one  that  minimizes 
the  maximum  value  his  opponent  can  obtain  against  it.  Con- 
versely, a best  strategy  for  the  maximizing  player  is  one  that 
maximizes  the  minimum  value  his  opponent  can  obtain  against 
it.  The  theory  of  mathematical  games  guarantees  that  such 
strategies  exist,  although  they  need  not  be  unique,  and  that 
the  nun-max  value  is  equal  to  the  max-min  value,  which  is 
called  the  value  of  the  game.  If  X*  and  Y*  are  used  to  denote 


solution  strategies  and  V*  the  value  of  the  game,  then 

V*  = min  max  V(X,i') 

X Y 

- max  min  V(X,Y) 

Y X 

= V(X*,Y*)  . 

Another  way  of  expressing  the  minimax  conditions  is 


V(X*,Y)  £ V*  for  all  Y 
V(X,Y*)  £ V*  for  all  X , 

which  is  to  say  that,  if  either  player  uses  a solution 
strategy,  the  other  player  cannot  find  a better  strategy 
than  his  own  solution  strategy. 

B.  Equivalent  Dual  Linear  Programs 

Charnes  (Reference  4)  has  shown  that  a constrained  game 
is  equivalent  to  dual  linear  programs  in  which  the  strategy 
vectors  are  augmented  by  components  that  are  Lagrangian 
multipliers.  In  our  case,  the  added  components  are  denoted 
by  superscript  zeros  or  subscript  zeros.  The  primal  program 
can  be  written: 


Maximize 


2 x°  D + l Xg  T 
„ n n _ o g 
n g ^ 


subject  to:  I l x?  (T  a?  ) = A 
i x g m 

Z x?  = 1 


Z x°  6g.  + x^  + Z x?  v? . < 0 


n nj 


i 13 
x?  > 0 
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where  all  indexes  are  to  be  read  as  greater  than  zero  unless 
explicitly  stated  as  zero.  Similarly,  the  dual  program  is: 


Minimize 
subject  to: 


Y a? 
m 


y° 
im  m 


Z A y°  + E T yg 
m 'm  g g -'o 


m 


Y.  7.  (T  6g.)  y?  * 
g j g n3  3 


* ♦ j v?j  y?  2 


t 


1 

0 


0 . 


Charnes  shows  that  solutions  of  the  dual  programs  exist  and 
that  the  value  of  the  equivalent  game  is 


V* 


Z x°  D 
n n n 


V 

g 


T 

g 


Y A y°  - Z T yg 
m m 7m  „ g o 
m 9 


From  this  result,  it  is  natural  to  interpret  - x°  and  - y° 
as  marginal  values  or  values  per  unit  of  the  respective 
resource  types,  and  to  interpret  - xg  and  - yg  as  intercept 
values  per  object  of  the  respective  object  types. 

C.  The  Q~Basis 

Balinski  and  Tucker  (Reference  5)  display  a scheme  for 
a Charnes-type  formulation.  We  suppose  such  a scheme  to 
exist  for  our  problem.  We  then  suppose  that  a solution  is 
known,  X*  and  Y*,  and  we  delete  from  the  scheme  all  those 
rows  and  columns  for  which  solution  components  x?  and  yg 
are  zero.  The  remaining  elements,  arranged  to  fit  the  block- 
diagonal  structure  of  the  problem,  are  illustrated  in 
Figures  3 and  4 as  what  we  call  a Q-basis.  It  is  a 
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representation  of  the  equations  satisfied  by  the  dual  solutions, 
considering  the  inequalities  as  side  conditions. 

In  Figure  3,  X and  Y are  the  augmented  solution  strategy 
vectors,  Q is  the  matrix  of  coefficients  of  the  basis 
equations,  and  A and  D are  the  right-hand  vectors.  In  matrix 
notation,  the  basis  equations  are  QY  = D and  XQ  = A. 

Figure  4 illustrates  the  partitioning  of  the  Q-basis  into 
blocks.  Superscripts  indicate  the  object  type  to  which  a 
block  pertains,  with  a superscript  zero  indicating  blocks 
associated  with  resources.  There  is  a band  of  resource- 
associated  blocks  across  the  top  of  X,  Q,  and  D and  a band 
down  the  left  side  of  Y,  Q,  and  A.  Down  the  main  diagonal 
of  Q is  a line  of  object-type  blocks.  The  rest  of  Q consists 
of  zero  elements. 

Figure  5 illustrates  the  elements  within  typical  resource 
blocks  and  object- type-g  blocks.  For  convenience,  the 
included  strategies  are  given  new  i and  j indexes  to  match 
their  positions  in  Q,  regardless  of  what  i and  j indexes 
they  might  have  had  originally.  The  elements  shown  as  zero 
are  always  zero  because  of  their  positions;  other  elements 
will  be  zero  only  if  one  of  their  variable  factors  is  zero. 

At  the  bottom  of  Figure  5 is  shown  a scheme  for  relating 
the  conditions  on  the  excluded  strategies  to  the  structure 
of  the  Q-basis.  If  the  vector  shown  is  denoted  by  I,  then 
the  matrix  condition  is  IY  i 0.  In  similar  fashion,  an 
excluded  maximizing  strategy  can  be  represented  by  a vector 
J and  the  exclusion  condition  by  XJ  < 0. 

The  various  blocks  of  the  Q-basis  must  be  conformable 
as  illustrated,  but  there  is  no  a priori  restriction  on 
their  size,  with  the  exception  that  every  Q^  block  must  be 
at  least  2x2,  since  at  least  one  elementary  strategy  for 
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each  player  is  required,  even  if  it  is  only  the  null  strategy. 
Otherwise,  Qg  may  have  any  dimensions  and  may  be  square  or 
rectangular,  and  if  rectangular  its  long  dimension  may  be  in 
either  direction.  In  general,  different  object  types  will 
have  Qg  blocks  of  different  size  and  shape. 

Likewise,  the  tot^l  Q matrix  itself  may  be  either  square 
or  rectangular,  but  it  is  nearly  always  square,  since  if 
Q were  to  be  rectangular,  one  player  would  have  more  variables 
in  his  augmented  strategy  than  conditions  for  them  to  satisfy. 
In  that  case  he  could  reasonably  adjust  to  a better  strategy 
restoring  the  balance  between  variables  and  conditions. 

Detailed  computations , involving  the  Q-basis  elements, 
show  that  the  value  of  the  game  is 

V*  = - AY* 

= - X*D  . 


Similarly,  the  minimax  conditions  can  be  verified  by 
using  the  side  conditions  on  excluded  strategies.  Further- 
more, if  Q is  square  and  non-singular,  then 

X*  = aq"1  , 

Y*  = Q_1D  , and 
v*  = - AQ-1D 


These  formulas  provide  a method  for  computing  the  solution 
if  the  Q-basis  is  known.  Section  III  will  describe  a method 
of  finding  the  basis. 


Before  leaving  the  discussion  of 

their  physical  feasibility  should  be 

ponents , x?  and  y?,  are  fractions  by 
t j 


solution  strategies, 
discussed.  The  corn- 
definition.  In  general. 
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the  products,  T xf  and  T y?,  will  still  be  fractions  and 
hence  not  physically  feasible.  Yet,  if  the  solution  strate- 
gies are  interpreted  as  mixed  strategies  in  a game-theoretic 
sense,  it  is  usually  possible  to  find  not  only  one,  but  many, 
equivalent  sets  of  physically  feasible  integral  strategies, 
with  an  associated  frequency  for  each  strategy  of  the  set. 

In  another  view,  the  fractions,  x?  for  example,  may  be 
considered  as  a frequency  distribution  on  each  object  of 
type  g,  with  the  important  proviso  that  the  distributions  on 
the  several  objects  of  type  g are  not  independent  of  each 
other  or  of  the  distributions  on  objects  of  other  types. 
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III.  PATH  METHOD 

The  path  method  of  solution  is  a form  of  parametric 
linear  programming  involving  variation  of  resources.  In  this 
section  we  will  describe  how,  as  resources  are  varied,  a Q- 
basis  can  be  changed  so  that  it  will  always  provide  a solution 
of  the  problem  for  the  current  numbers  of  resources. 

A.  Space  of  Resources 

It  is  helpful  to  think  of  any  set  of  resources 

(Aj , * • • , am,*  Dj  , • ♦ * , nN) 

as  defining  a point  in  an  (M  + N) -dimensional  space  of 
resources.  Associated  with  each  point  of  this  space  is  a 
distinct  problem,  or  matrix  game,  which  has  at  least  one  pair 
of  solution  strategies,  X*  and  Y*,  and  a unique  value  of  the 
game , V* . 

The  values,  V*,  may  be  thought  of  as  plotted  in  (M  + N + 1) 
dimensions  so  as  to  form  a continuous  solution  surface  over 
the  space  of  resources.  At  nearly  every  point  of  this  surface 
there  is  a unique  tangent  hyperplane,  whose  slopes  in  the 
coordinate  directions  are  the  same  as  the  marginal  values 
defined  in  Section  II,  i.e., 

, o o o o, 

<-yi,  •••»  -yM;  -xlt  •••,  -xN) 

For  each  point  of  the  resource  space  there  is  at  least 
one  valid  Q-basis  from  which  solution  strategies  and  the 
value  of  the  game  may  be  computed.  There  may  be  more  than 
one  such  valid  Q-basis,  each  yielding  different  solution 
strategies,  but  there  is  only  one  value  of  the  game  at  any 
point. 
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For  nearly  every  point  of  the  resource  space ,.  a valid 
Q matrix  is  square  and  non-singular.  Such  a basis  is  usually 
valid  for  a continuous  set  of  points  that  we  will  call  a 
region  of  the  resource  space.  Throughout  a region,  the 
strategies  are  uniquely  defined  by  the  basis. 

The  interior  of  a regio::  is  generally  characterized  by 
proper  inequalities: 

. x?  > 0 for  included  i 

• y?  > 0 for  included  j 

, IY  > 0 for  excluded  i 

XJ  < 0 for  excluded  j 

The  boundary  of  a region  is  characterized  by  one  or  more  of 
these  quantities  becoming  equal  to  zero. 

fl . The  Path 

The  continuous  variation  of  resources  from  some  initial  set 
of  values  in  some  prescribed  fashion  to  some  terminal  set  of 
values  generates  a path  through  the  resource  space.  The  path 
method  enables  us  to  find  solutions  all  along  a fairly  arbi- 
trary path  if  we  have  a solution  basis  at  the  initial  point. 

There  is  one  point  in  every  type  of  problem  where  the 
solution  basis  is  always  known.  That  point  is  the  origin, 
where  all  resources  are  zero  and  both  players  use  null  strate- 
gies. Hence,  a path  can  always  start  at  the  origin,  as  it 
does  in  the  computer  programs  of  Sections  IV  and  V.  A path 
may,  however,  start  at  any  other  point  where  the  Q-basis  is 
known. 

The  simplest  kind  of  path  is  generated  by  varying  only 
one  resource  type  at  a time,  all  the  others  being  held  con- 
stant. We  will  call  this  a rectangular  path,  since  it  con- 
sists of  a series  of  straight-line  segments  at  right  angles 
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to  each  other  in  the  resource  space.  This  kind  of  path  is 
used  in  the  programs  of  Sections  IV  and  V.  In  this  method, 
the  dimensions  of  the  problem  start  at  (M,N)  - (0,0)  and 
increase  as  resource  types  are  introduced.  Once  a type  has 
been  introduced,  the  number  of  resources  of  that  type  may  be 
increased  or  decreased  at  will. 


A more  general  kind  of  path  is  one  consisting  of  straight- 
line  segments  generated  by  simultaneously  varying  any  or  all 
of  one  player's  resources  linearly  in  terms  of  some  parameter. 
A segment  of  such  a path  can  be  represented  by  a set  of 
equations  of  the  form 


Am  = A,«  + A,„,h'0ih-h  / 
m om  lm 

expressing  the  condition  that  the  number  of  type  m resources 
varies  from  an  initial  value  of  A.  when  the  parameter  h = 0 

OKI 

to  a terminal  value  of  + A h.  In  similar  fashion,  a 

o m x m 

variation  of  maximizing  resources  can  be  expressed  as 


D = D + D h . 
n no  m 

The  parameter  h may  be  arbitrary  or  it  may  have  some  real 
meaning,  e.g.,  a vehicle  that  carries  fixed  numbers  of  each 
resource  type,  or  a budget  allocated  in  fixed  proportions  to 
the  different  resource  types.  Either  minimizing  or  maximizing 
resources  may  be  varied  on  a path  segment,  but  we  have  no 
convenient  method  of  varying  both  simultaneously. 


C.  Location  of  a Regional  Boundary 

If  a path  segment  lies  entirely  in  a single  region,  the 
same  Q-basis  provides  a solution  at  every  point  of  the  segment. 
If  not,  the  first  problem  is  to  find  the  value  of  h for  which 
the  path  meets  a regional  boundary.  This  critical  value  of  h 
defines  a critical  point  on  the  path  where  some  changes  must 
be  made  in  the  Q-basis  before  proceeding. 
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To  do  this  for  a minimizing  path,  for  example,  we  first 
add  a row  to  the  A vector  of  Figure  4 to  reflect  the  variable 
nature  of  the  resources,  i.e.,  A becomes  the  two-row  matrix 


Aei'  Aom'TJ'0'  Tg'°'  ° 

AM,  ••*,  A^,0  #0,  ••*,  0 ,0,  •••,  0 


assumed  to  be  multiplied  by  H = (1,  h) . The  matrix  equation 
of  condition  then  becomes 


XQ  s HA  . 

The  vector  solution 

X = HAQ~ 1 

is  then  equivalent  to  a set  of  scalar  equations  of  the  form 

x = x0  + xth 

One  way  the  basis  can  fail  is  by  some  x|  for  included  i 
becominq  negative.  If  xt  > 0,  x cannot  become  negative  as 
h increases  and  we  pass  on  to  the  next  x.  If  x;  < 0,  we 
solve  the  equation 

h = -x0/x, 

to  find  the  value  of  h for  which  x = 0.  The  least  of  these 
values  of  h is  a candidate  for  the  critical  value. 

The  second  way  the  basis  can  fail  is  for  XJ  to  become 
greater  than  0 for  some  excluded  j . The  scalar  product  XJ 
is  of  the  form 

xe  + Xj  h . 

Applying  a two-part  test  for  every  excluded  j,  we  pass  on 
to  the  next  if  xt  < 0,  since  then  x cannot  become  positive. 
If  Xj  > 0 , we  solve  the  equation 

h - -Xj/X! 
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to  find  the  value  of  h for  which  XJ  » 0.  The  least  of  these 
values  is  a second  candidate  for  the  critical  value. 

The  lesser  of  the  two  candidates  is  the  critical  value 
of  h,  which  by  substitution  defines  the  critical  point  on 
the  path  and  the  solution  strategy  for  the  minimizing  player 
at  that  point. 

If  the  critical  value  comes  from  x?  = 0,  then  the 
corresponding  row  of  the  Q-basis  must  be  deleted;  if  it 
comes  from  XJ  = 0,  then  a corresponding  column  must  be  added 
to  the  Q-basis.  In  either  case  the  resulting  Q-basis  will 
have  one  more  column  than  it  has  rows. 

In  similar  fashion,  with  appropriate  formulations  and 
test  criteria,  we  can  find  the  critical  value  along  a path 
in  the  space  of  resources  of  the  maximizing  player;  in  that 
case,  the  resultant  change  in  the  Q-basis  is  deletion  of 
a column  or  addition  of  a row,  and  the  Q-basis  will  have  one 
more  row  than  it  has  columns. 

D.  Finding  the  Q-Basis  for  the  Next  Region 

The  rectangular  Q-basis  resulting  from  the  location  of  a 
boundary  is  a valid  basis  on  the  boundary.  Taking  the 
example  of  a minimizing  path,  there  is  one  more  column  than 
row,  and  hence  one  more  condition  on  the  minimizing  player 
than  he  has  variables.  Although  the  extra  condition  is 
redundant  and  is  satisfied  by  the  minimizing  player's 
solution  strategy  at  the  critical  point,  it  does  act  as  an 
added  constraint. 

On  the  other  hand,  the  maximizing  player  here  has  one 
more  variable  than  conditions  to  be  satisfied.  Because  of 
this,  there  is  an  infinite  number  of  solutions  for  the 
maximizing  player.  In  fact,  any  linear  combination  of  his 
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solutions  on  the  two  sides  of  the  boundary  is  a solution 
for  him  on  the  boundary. 

it  is  the  degree  of  freedom  of  the  rectangular  Q-basis 
that  lets  the  maximizing  player  shift  his  strategy  on  the 
boundary  to  that  strategy  appropriate  in  the  next  region. 
We  determine  the  new  strategy  by  means  of  a new  parameter, 
the  marginal  value  per  unit  of  the  minimizing  resource 
parameter  and  denote  the  new  parameter  by  e,  defined  as 


e 


-A  V° 
A1  Iyj 


AlMyM 


We  represent  this  condition  by  temporarily  adding  in 
the  last  row  of  the  Q matrix  the  vector 


(A 


1 1 


0) 


adding  to  the  D matrix  a second  column  consisting  of  all  Os 
except  for  -1  in  the  last  row,  and  assuming  D to  be  multi- 
plied by  the  vector  E - ^ | . The  matrix  equation  of  con- 
dition becomes 


QY  = DE 


The  vector  solution 


Y - Q-1  DE 

is  equivalent  to  the  set  of  scalar  equations  of  the  form 

y - y0  + yie  . 


The  maximizing  player  is  interested  in  driving  e to  a 
value  as  high  as  possible  algebraically,  that  is  in  maxi- 
mizing survival  in  the  direction  of  the  path.  But  he  is 
limited  by  the  conditions,  y?  > 0 for  every  included  j,  and 
IY  > 0 for  e^ery  excluded  i.  The  critical  value  of  e is  the 
lowest  value  meeting  one  of  these  conditions. 
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Applying  a two-part  test  to  y,  we  pass  on  if  yt  £ 0,j 
since  y cannot  then  become  negative  as  e increases.  If 
yj  < 0,  we  solve  the  equation 

e - -y./y, 

to  find  the  value  of  e for  which  x - 0.  The  least  of  these 
values  of  e is  a candidate  for  the  critical  value. 


The  second  test  is  on  the  condition  IY  > 0.  This  product 
is  of  the  form 


y0  + y 

We  pass  on  to  the  next  i if  y,  i 0,  since  y cannot  then 
become  negative  as  e increases.  If  yj  < 0,  we  solve 

e = -v6/yl 

and  select  the  least  of  these  values  as  a second  candidate 
for  a critical  value. 

The  lesser  of  the  two  candidates  is  the  critical  value 
of  e,  which  by  substitution  defines  the  maximizing  strategy 
along  the  path  into  the  next  region. 

The  rectangular  Q-basis  is  then  changed  by  deleting  the 
column  or  adding  the  row  associated  with  critical  e.  This 
change  restores  Q to  a square,  non-singular  matrix,  the 
solution  basis  for  the  next  region. 

In  similar  fashion,  with  appropriate  formulations  and 
tests,  we  find  the  Q-basis  for  the  next  region  along  a 
maximizing  path. 

E.  Terminating  a Path  Segment 

The  two  processes  described  in  Sections  III.C.  and 
III .0.  repeat  in  alternation  until  the  end  of  a path  segment 
is  signaled  by  critical  h Z h.  At  this  interior  point  of  a 
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region,  the  Q-basis  is  not  changed,  the  terminal  strategy 
can  be  determined  by  substitution  of  h,  the  terminal  coordi- 
nates are  entered  in  the  first  row  of  A or  first  column  of 
D,  as  the  case  may  be,  and  either  a new  path  segment  is 
begun  or  the  problem  is  ended. 

F.  Beginning  a Path  Segment 

The  previous  exposition  was  given  with  M + N resource 
types  being  allocated.  If  the  new  path  segment  specifies 
a variation  of  no  more  than  these  resource  types,  then  the 
terminal  point  of  the  previous  segment  is  an  interior  point 
of  the  space  and  we  begin  the  new  path  by  finding  a boundary 
as  in  Section  1XI.C.  above. 

If,  however,  a new  type  of  resource  is  introduced,  the 
terminal  point  of  the  previous  segment  automatically  becomes 
a boundary  point  of  the  new  (M  + N + 1) -dimensional  space  of 
resources.  The  Q-basis  must  be  expanded  by  opening  a row 
or  column  of  zeros  at  the  appropriate  place.  The  new  path 
must  then  be  begun  by  using  the  process  described  in  Section 
III.D. 

G.  Beginning  a Problem 

Since  the  method  will  not  permit  introducing  more  than 
one  type  of  resource  at  a time,  it  is  usual  to  begin  a 
problem  at  the  origin,  where  both  minimizing  and  maximizing 
players  are  constrained  to  a null  strategy  and  where  the 
Q-basis  has  2G  rows  and  columns  as  shown  in  Figure  6. 

The  problem  is  then  begun  by  introducing  one  resource 
type  in  any  desired  numbers,  followed  by  the  introduction 
of  other  resource  types  until  all  M + N have  been  intro- 
duced . 
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FIGURE  6 ; 

Q-BASIS  AT  THE  ORIGIN  j 

The  order  of  introduction  may  make  some  difference  in  the  j 

computational  time  but  makes  no  difference  in  the  terminal  ' 

value  of  the  game.  However,  there  are  certain  simple  j 

, . s 

restrictions  that  must  be  observed.  For  example,  the  marginal  1 

value  must  not  be  equal  to  zero.  This  would  be  the  case  in  ’ 

the  weapon  allocation  problem  if  we  tried  to  introduce  defense 

before  attack,  since  in  that  case  all  targets  would  survive 

j 

regardless  of  whatever  allocation  the  defense  might  make. 
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IV.  PATHS?  COMPUTER  PROGRAM 


The  PATHS 7 computer  program  is  explained  in  this  section. 
It  is  specifically  designed  to  solve  the  multi -dimensional 
two-sided  weapon  allocation  game  for  point  targets.  However, 
it  can  be  modified  with  little  trouble  to  solve  a variety 
of  resource  allocation  problems.  In  fact,  ease  of  modifica- 
tions has  been  one  of  the  major  criteria  influencing  program 
design. 

The  program  is  written  4 a BASIC  language  for  use  on  an 
IBM  360/65  computer  in  a time-sharing  mode  with  interaction 
of  computer  and  operator.  The  computer-system  constraints 
that  have  been  binding  at  one  time  or  another  during  the 
evolutionary  development  of  the  program  are: 

-1  A limit  of  800  statement  lines 

-2  A limit  of  80  FOR  loops 

-3  A limit  of  29  numeric  arrays 

-4  Limited  storage  space  for  the  array  elements  of 
the  problem. 

The  first  three  of  those  constraints  have  been  overcome 
by  such  devices  as  using  the  same  subroutines  to  serve  both 
the  minimizing  and  the  maximizing  player  and  using  the  same 
arrays  to  store  similar  numbers  associated  with  both,  thus 
taking  advantage  of  the  structural  symmetry  of  the  problem 
and  effectively  transposing  large  matrices  back  and  forth 
by  just  changing  a few  indices. 

The  effect  of  the  fourth  constraint  has  been  minimized 
by  using  a subroutine  to  compute  the  value  function  as  needed 
instead  of  precomputing  and  storing  it.  Also,  some  matrices 
have  been  reduced  in  size  by  packing  the  significant  elements. 

Preceding  page  blank 
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Computer  processing  time  has  been  reduced  by  designing 
a three-stage  solution  process  with  recursions.  Further 
improvements  in  processing  time  have  resulted  from  the  use 
of  single  indices  instead  of  double  for  nearly  all  arrays. 

The  price  paid  for  the  increase  in  capacity  and  decrease 
in  running  time  has  been  some  rather  complicated  indexing. 
However,  even  that  is  not  completely  without  value,  since 
the  indexing  is  an  aid  to  flexibility. 

A.  Program  Structure 

The  program  is  composed  of  an  initialization  section, 
a control  section,  and  a collection  of  subroutines.  The 
heart  of  the  program  is  the  control  section,  which  calls 
the  main  subroutines  as  needed. 

Blocks  of  statement  lines  are  allotted  as  follows: 


(1-999) 

(1000-1999) 

(2000-2999) 

(3000-3999) 

(4000-4999) 

(5000-5999) 

(6000-6999) 

(7000-8159) 

(8160-8999) 

(9000-9999) 


Initialization  section 

Control  section 

Test  subroutine 

Auxiliary  subroutines 

Value  subroutine 

Print  subroutine 

Addition  subroutine 

Deletion  subroutine 

Auxiliary  subroutines 

Strategy  subroutine. 

The  remainder  of  this  section  contains  a complete  listing  of 
the  program  statements  in  numerical  order,  with  explanatory 
comments  following  each  small  group  of  statements. 
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Initialisation  Section 


This  section  sets  the  dimensions  of  the  problem, 
initializes  various  indices,  inputs  data,  and  makes  pre- 
liminary computations. 


L 
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These  lines  provide  an  example  of  the  data  that  an  operator 
must  type  lor  any  problem  he  proposes  to  run.  Other  examples 
appear  in  Section  VI,  Cases  5 to  8. 

The  first  four  numbers  define  the  dimensions  of  the 
problem.  In  this  example,  there  are  2 target  groups  (object 
types),  of  which  1 is  defended,  2 attack-weapon  types,  and 
2 defense- veapon  types.  The  program  later  reads  these 
numbers  as  Gl,  G4,  Al,  and  Dl,  respectively.  By  convention, 
the  defended  target  groups  are  the  first  G4  groups. 

The  next  2*G1  numbers  define  target  data,  that  is  the 
number  of  targets  and  value  per  target  for  each  group. 

The  next  G1*A1  numbers  define  attack-weapon  data  as  a 
matrix  of  single-shot  kill  probabilities,  in  order  of  all 
weapon  types  against  the  first  target  group,  then  all  weapon 
types  against  the  second  target  group,  etc.  The  weapon  types 
must  be  listed  in  the  same  order  for  each  target  group,  the 
preferred  order  being  from  least  effective  weapon  to  most 
effective  weapon,  as  explained  in  Section  IV. F. 

The  last  Dl  numbers  define  defense—  weapon  data  as  single- 
shot probabilities  of  intercept,  it  being  assumed  that  these 
are  the  same  against  every  type  of  attack  weapon.  The 
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preferred  order  of  listing  is  from  roost  effective  defense 
weapon  to  least  effective,  as  explained  in  Section  IV. F. 
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These  lines  reserve  storage  space  in  core  for  23  lists  and 
arrays,  those  with  related  dimensions  appearing  on  the  same 
line.  The  reservations  are  adequate  for  most  problems 
involving  no  more  than  100  target  groups  (Gl)  and  no  more 
than  15  weapon  types  in  all  (A1  + Dl)  . 

All  of  the  matrices  will  now  be  discussed,  in  a con- 
venient order.  The  form  (T)  will  denote  a matrix  itself, 
and  the  forms  T(2) , T (G) , $ (A7  + 1),  etc.,  will  denote 
particular  elements  of  a matrix. 

(T)  stores  the  number  of  targets  in  each  group.  It  is 
customarily  addressed  by  the  simple  variable  G,  i.e.,  T(G). 

(V)  stores  the  value  per  target  in  each  group  and  is 
usually  addressed  as  V(G). 

(P)  stores  the  probabilities  of  kill.  Conceptually,  (P) 
is  a two-dimensional  array,  but  it  is  treated  by  the  program 
as  a one-dimensional  list.  The  manner  of  address  is  discussed 
in  Section  IV. F.  in  connection  with  the  value  function, 

SUB  4000. 

(D)  stores  the  complements  of  the  probabilities  of  inter- 
cept. The  manner  of  address  is  discussed  in  connection  with 
the  value  function. 

($)  stores  tho  force  levels  attained  at  the  end  of  any 
path  segment.  Elements  $(1)  to  $ (Al)  represent  levels  of 
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attack  weapons/  and  $ (A1  + 1)  to  $ (Al  + Dl)  levels  of 
defense  weapons. 

(A)  is  the  largest  matrix  in  the  program.  Primarily, 
it  stores  those  elementary  strategies  of  both  attack  and 
defense  that  are  in  the  basis  at  any  time.  It  also  stores 
some  other  numbers  related  to  those  elementary  strategies. 
Conceptually,  it  has  an  internal  structure  that  is  illustrated 
by  the  storage  map  of  Figure  7.  The  matrix  is  partitioned 
into  three  main  cell  regions:  test  (1-60),  attack  (61-3000), 

and  defense  (3001-4000) . The  test  region  will  be  discussed 
in  connection  with  the  testing  process,  SUB  2000.  The 
attack  region  is  subdivided  into  rows,  each  row  having 
Al  + 1 elements.  The  example  showing  6 elements  is  based 
on  Al  = 5.  These  rows  are  grouped.  The  first  group, 
consisting  of  a single  row,  stores  the  identifying  numbers 
of  the  attack  weapons  currently  in  play,  e.g. , 2,  3,  and  5. 

The  second  group  of  G1  rows  stores  test  control  data  derived 
from  the  elementary  strategies.  The  remaining  rows  store 
the  elementary  strategies  for  each  target  group  in  vector 
form.  For  example,  the  second  elementary  strategy  on  target 
group  3 has  0 type- 2 weapons,  2 type- 3 weapons,  and  1 type-5 
weapon,  for  a total  of  3 weapons.  The  maximum  entries  that 
have  occurred  in  any  elementary  strategy  row  of  target  group  3 
during  the  course  of  the  run  are  stored  in  the  third  row  of 
the  test  control  group,  e.g.,  2,  2,  2,  and  4.  The  defense 
region  is  structured  and  used  like  the  attack  region,  but 
the  rows  have  Dl  + 1 elements,  three  in  the  example.  The 
address  of  an  element  of  (A)  is  usually  compounded  from  its 
place  in  its  row  and  the  address  of  the  terminal  element  of 
the  preceding  row,  e.g.,  A(120  + 4)  = 3. 
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(X)  and  (Y)  store  computed  solutions  of  the  parametric 
equations  of  the  Q~basis.  (X)  stores  the  first  column  of  a 
solution  and  (Y)  the  second  column,  i.e.,  the  column  which 
is  multiplied  by  the  parameter.  They  are  used  together, 
sometimes  for  an  attack  solution,  sometimes  for  a defense 
solution.  Whichever  the  case,  the  elements  of  (X)  and  (Y) 
are  arranged  in  the  appropriate  sequence  established  in 
Section  II  and  illustrated  in  Figure  5s  marginal  values  for 
opposition  weapons,  intercept  value  on  group  1,  elementary 
strategies  on  group  1,  intercept  value  on  group  G1 , 

elementary  strategies  on  group  Gl. 

Before  discussing  the  other  matrices  of  the  program,  it 
is  desirable  to  describe  in  general  terms  a three-stage 
solution  process  that  saves  storage  space  and  computer 
running  time.  Instead  of  storing  and  inverting  a Q-matrix 
whose  dimensions  would  exceed  200  x 200  if  there  were  100 
target  groups,  we  allow  the  Q-basis  of  Sections  II  and  III 
to  exist  only  as  a mental  concept  of  the  equations  of  the 
problem.  These  equations  are  then  solved  by  two  stages  of 
Gaussian  elimination,  a third  stage  in  which  a small  matrix 
is  inverted,  and  a back  solution  process  which  computes  (X) 
and  (Y)  as  full-size  solutions  of  the  conceptual  Q-basis. 
Intermediate  results  are  stored  in  such  form  that  they  can 
be  modified  recursively,  without  having  to  repeat  the  entire 
three-stage  process  eve  y time  a line  is  deleted  from  or 
added  to  the  conceptual  Q-basis.  In  the  program,  the 
matrices  (Q) , (R) , (S) , (U) , (W) , (Z),  (C) , and  (F)  are 
primarily  involved  in  the  process.  These  will  now  be 
discussed . 

(R)  stores  the  results  of  the  first-stage  Gaussian 
elimination,  which  reduces  the  Q-basis  by  2 rows  and  2 columns 
for  each  target  group.  The  process  is  precisely  defined.  In 


35 


the  example  of  Figure  5,  the  upper  left  corner  of  the  central 
block  hag  this  configuration— 

0 © 

© V'.  • 

Pivots  on  the  circled  elements  eliminate  two  rows  and 
columns.  Because  of  the  choice  of  pivot  elements r the 
computations  involve  only  subtractions,  and  vary  little 
information  is  lost  by  physically  eliminating  the  pivot 
rows  and  columns.  In  fact,  since  the  elementary  strategies 
are  preserved  in  (A)  and  the  number  of  targets  in  (T) , we 
need  only  provide  for  saving  vM  elsewhere.  (R)  must  be 
large  enough  to  hold  the  elements  not  physically  eliminated. 

Its  size  varies  during  a run.  At  the  beginning,  or  at  the 
origin,  (R)  is  zero  since  all  rows  and  columns  are  eliminated 
(see  Figure  6).  Then,  us  rows  and  columns  are  added  to  and 
deleted  from  the  conceptual  Q~basis,  (R)  is  modified  to 
reflect  the  changes.  Suitable  recursion  algorithms  are  in 
the  program  for  that  purpose. 

(Q)  stores  the  values  v,j  for  each  target  group.  These 
"base”  values  would  otherwise  be  lost  after  the  first-stage 
elimination. 

(S)  stores  the  results  of  the  second-stage  Gaussian 
elimination,  but  preserves  the  physical  configuration  of 
(R) . In  concept,  both  (R)  and  (S)  are  rectangular,  as 
illustrated  by  the  pattern  of  Figure  8.  However,  the 
elements  are  stored  according  to  the  numbering  scheme  in  the 
figure.  These  matrices  are  augmented  to  include  elements 
from  the  right-hand  sides  of  the  equations  of  the  problem. 

So  the  first  column  and  first  row  of  Figure  8 contain  elements 
from  vectors  labeled  D and  A in  Figure  3.  The  remaining 
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columns  and  rows  correspond  to  the  none limina ted  columns 
and  rows  of  the  Q-basis.  In  this  example,  the  columns  repre- 
sent 3 attack  weapon  types  and  5 excess  elementary  defense 
strategies,  and  the  rows  represent  2 defense  weapon  types 
and  6 excess  elementary  attack  strategies.  The  lower  right 
region  of  Figure  8,  which  we  will  call  the  pivot  region, 
contains  elements  left  over  after  the  first-stage  elimination 
of  two  rows  and  columns  from  each  of  the  blocks  on  the 
diagonal  of  the  Q-basis.  In  this  example,  supposing  there 
are  3 target  groups  in  the  problem,  the  gap  in  the  pivot 
region  implies  that  the  original  Q-basis  block  for  group  2 
has  3 rows  but  only  2 columns,  and  hence  is  not  explicitly 
represented  in  (R)  and  (S) . The  blank  elements  in  the  pivot 
region  are  understood  to  be  zeros,  which  need  not  be  stored. 
The  second-stage  Gaussian  elimination  may  be  described  as 
a partial  inversion  in  place,  since  some  of  the  procedures 
are  adapted  from  the  MATINV  subroutine  (Reference  6)  for  a 
total  inversion  in  place.  In  concept,  (R)  is  the  given 
matrix  and  (S)  is  the  result  of  the  partial  inversion.  As 
many  as  possible  pivot  elements,  such  as  those  identified  in 
Figure  S,  are  selected  from  the  pivot  region  only,  so  that 
each  pivot  affects  only  one  of  the  target  groups.  There  is 
no  shifting  of  rows  or  columns,  but  the  location  of  pivots 
is  recorded.  The  pivot  rows  and  columns  are  actually  used 
to  store  the  results  of  the  partial  inversion,  but  it  is 
understood  that  a pivot  row  or  column  has  an  alternative 
aspect,  characteristic  of  a Gaussian  elimination,  in  which 
the  pivot  element  is  1 and  all  the  other  elements  of  the  row 
or  column  are  0.  The  proper  aspect  is  used  at  any  point  of 
the  program.  Finally,  as  lines  are  changed  in  the  conceptual 
Q-basis,  recursion  algorithms  make  the  induced  changes  in  (S) . 
At  the  end  of  the  second  stage,  the  number  of  equations 
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remaining  to  be  solved  is  somewhat  unpredictable,  but  in 
theory  can  be  bounded:  Max  {Al,Dl}  < number  < A1  + Dl. 

In  the  example  of  Figure  8,  four  equations  remain  to  be 
solved  out  of  an  original  fourteen  in  the  conceptual  Q-basis. 

(0),  (W) , and  (Z)  are  used  in  the  third  stage,  which 
solves  the  remaining  equations  by  matrix  inversion.  On  each 
occasion,  they  are  redimensioned  by  the  program  to  precisely 
the  size  needed,  e.g.,  in  the  continuation  of  Figure  8 they 
are  redimensioned  as  4 x 4 matrices.  (U)  is  then  filled 
with  the  elements  of  (S)  that  are  in  neither  a pivot  row 
nor  a pivot  column,  as  in  Figure  9. 
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FIGURE  9 

FORMATION  OF  MATRIX  <U)  FROM  ELEMENTS  OF  MATRIX  (S) 


(W)  stores  the  inverse  of  (U) . (Z)  stores  either  a duplicate 

or  the  transpose  of  (W) , whichever  is  needed  at  the  time. 
These  three  matrices  are  the  only  ones  with  two  subscripts. 

(C)  and  (F)  store  the  two-column  solution  of  the  third- 
stage  equations.  Their  elements  become  elements  of  (X)  and 
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(V) , respectively,  and  are  used  in  the  back  solution  to 
complete  these  vectors. 

(G)  is  the  first  of  seven  index  matrices  that  are  needed 
to  help  locate  numbers  in  the  primary  matrices  already  dis- 
cussed. (G)  stores  the  locations  in  (S)  of  the  initial 
elements  of  blocks  of  numbers  pertaining  to  the  different 
target  groups.  (G)  consists  of  three  regions,  each  having 
G1  + l elements.  The  first  region  contains  the  locations 
of  blocks  in  the  pivot  region  of  (S) . For  the  example  of 
Figure  8,  G(l)  = 1,  G(2)  = 7 since  this  is  where  the  block 
for  the  second  target  group  would  appear  if  it  existed, 

G ( 3)  = 7,  and  G(4)  = 13,  in  effect  defining  an  upper  bound 
on  the  pivot  region.  Similarly,  the  next  four  elements 
locate  blocks  in  the  lower  left-hand  region  of  (S) , i.e., 

G (5)  = 401,  G (6)  = 409,  G{7)  « 413,  G(8)  = 425.  The  last 
four  elements  locate  blocks  in  the  upper  right  region  of 
(S),  i.e.,  G (9)  * 751,  G(10)  = 760,  G(ll)  = 760,  G(12)  = 766. 
In  a problem  with  only  three  target  groups,  the  remaining 
elements  of  (G)  would  never  be  used. 

(N)  stores  the  number  of  rows  and  columns  of  each  block 
in  the  pivot  region  of  (S) . (N)  consists  of  two  regions  of 

G1  + 1 elements  each,  the  first  region  for  rows,  the  second 
for  columns.  So,  when  (5)  is  as  shown  in  Figure  8,  this 
index  will  be: 

(N)  - (2,1,3,0?3,0,2,0;0,« ••)  . 

Actually,  the  G1  + first  element  of  each  region  is  super- 
fluous, but  is  retained  for  convenience  in  indexing  (N) 
itself. 

(M)  stores  the  total  number  of  rows  and  columns  of  the 
pivot  region  preceding  each  block  of  the  pivot  region  of  (S) . 
(M)  is  arranged  the  same  way  as  (N) . For  our  standard 
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example,  we  have: 

(M)  s (0, 2 , 3, 6 ; 0 , 3 , 3, 5 } 0 , • • • ) 

In  this  case,  the  first  element  of  each  region  might  be 
considered  superfluous. 

(0)  stores  the  locations  of  the  actual  pivot  elements 
in  (S) , identifying  them  by  rows  and  columns  rather  than  by 
a linear  storage  number.  (0)  is  divided  into  two  regions  of 
equal  size,  in  this  case  having  30  elements  each.  The  first 
region  of  0 contains  an  entry  for  each  row  of  the  pivot  region 
of  (S) . If  the  row  has  no  pivot  element  in  it,  the  entry  is 
0.  If  the  row  has  a pivot  element  in  some  column,  the  entry 
is  the  number  of  that  column  within  the  block  of  columns  for 
the  particular  target  group  involved.  For  the  example  of 
Figure  8,  we  have  0(1)  to  0(6)  given  as  (1 , 2 , 0 , 2 ,0 , 1) . The 
second  region  contains  similar  entries  that  identify  a row 
for  any  pivot  in  a column,  so  we  have  0(31)  to  0{35)  given 

as  (1,2, 0,3,1) . 

(1)  is  the  central-index  matrix.  It  stores  a variety 
of  constants  and  variables  that  will  be  individually  dis- 
cussed in  connection  with  lines  140-215.  Its  conceptual 
structure,  shown  in  Figure  .10,  reflects  the  paired  nature 

of  most  of  its  elements.  It  has  a central  spine  of  elements 
identified  (except  for  1)  by  numbers  of  the  form  3n,  and 
two  wings,  of  the  forms  3n-l  and  3n+l.  The  elements  of  one 
wing  are  associated  with  rows  of  the  problem  and  those  of 
the  other  wing  with  columns.  An  interchange  of  wings  has 
the  effect  of  transposing  the  0-basis  of  the  problem. 

(K)  and  (L)  provide  the  mechanism  for  the  transposition, 
one  wing  of  (I)  being  read  into  (K)  and  the  other  into  (L) 
by  SUB  8690.  A flag  Si,  whose  value  is  set  at  +1  or  -1, 
determines  which  wing  is  read  into  which  matrix.  Most  of 
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the  computational  algorithms  of  the  program  are  written  in 
terms  of  (K)  and  (L)  and  are  equally  valid  for  row  operations 
or  column  operations,  that  is  for  attack  strategies  or 
defense  strategies. 
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These  lines  read  and  print  the  dimensions  of  the  problem, 
i.e.,  the  data  of  line  50. 
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This  sequence  computes  and  initializes  various  indices,  all 
others  being  automatically  initialized  at  zero  by  the  BASIC 
system.  The  paired  structure  of  (I)  is  illustrated  in  the 
following  discussion. 

K (1) , L (1 ) : I (2)  is  the  number  of  types  of  defense  weapons 

already  brought  into  play,  sometimes  denoted  by  D2.  1(4) 

is  the  number  of  types  of  attack  weapons  already  brought  into 
play,  sometimes  denoted  by  A2.  Both  1(2)  and  1(4)  are 
variables,  automatically  initialized  at  0 and  stepped  as  new 
weapon  types  are  added.  K(l)  and  L(l)  are  used  principally 
as  upper  limits  for  loops  and  as  locators  in  (X)  and  (Y) . 
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K (2) , L(2)  : 1(5)  and  1(7)  are  constants,  Dl  and  Al, 

respectively.  They  are  not  used  in  this  version  of  the 
program. 

K (3) , L(3) : 1(8)  and  1(10)  are  constants,  0 and  Al, 

respectively.  They  are  used  as  locators  for  the  two 
regions  of  the  ($)  matrix. 

K (4) , L (4) ; Not  used  in  this  version  of  the  program. 

K (5) , L (5)  ; 1(14)  and  1(16)  are  constants,  0 and  30, 

respectively.  They  are  used  as  locators  for  the  two 
regions  of  the  (0)  matrix  or  pivot  index.  In  case  (0) 
is  redimensioned,  for  example  to  80,  then  1(16)  should  be 
set  at  1/2  Dim (0)  * 40. 

K(6) , L (6) : 1(17)  and  1(19)  are  constants,  set  at  computed 

values  G9  = Gl  + 1 and  G8  = 2*G9,  respectively.  They  are 
used  as  locators  for  the  two  upper  regions  of  the  (G)  matrix. 

K(7) , L(7) : 1(20)  and  1(22)  are  constants,  0 and  G9, 

respectively,  serving  as  locators  for  the  two  regions  of 
the  (M)  and  (N)  matrices. 

K ( 8) , L ( 8) ; I (23)  and  1(25)  are  constants,  used  as  locators 
for  the  attack  and  defense  regions,  respectively,  of  (A). 

In  this  version  of  the  program  they  are  set  arbitrarily  at 
60  and  3000.  The  apportionment  of  space  between  attack  and 
defense  can  be  modified  by  changing  1(25). 

K ( 9 ) , L ( 9 ) : 1(26)  and  1(28)  are  constants,  used  to  specify 

the  number  of  elements  of  (A)  needed  to  record  the  useful 
features  of  an  elementary  strategy  for  the  attack  or  defense, 
respectively.  In  this  version  of  the  program  the  number  of 
elements  is  one  plus  the  number  of  weapon  types,  so  that 
the  number  of  weapons  of  each  type  plus  the  sum  of  the 
weapons  of  all  types  can  be  stored.  The  effect  is  to 
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structure  the  attack  region  of  (A)  into  rows  of  size  1(26) 
and  the  defense  region  into  rows  of  size  1(28). 

K ( 10) , L(10) ; 1(29)  and  1(31)  are  variables,  used  to  specify 

the  number  of  elements  currently  in  use  in  each  of  the  rows 
just  described.  Both  are  initialized  at  the  same  value  as 
1(27).  Thereafter,  1(29)  is  kept  equal  to  D2  + 1(27)  and 
1(31)  to  A2  + 1(27) . 

K (11) , L ( 11 ) ; 1(32)  and  1(34)  are  constants,  locating  the 

initial  values  of  the  a and  6 regions  of  (5)  respectively  and 
are  set  at  401  and  751  in  this  version  of  the  program.  They 
may  be  changed  if  a different  apportionment  of  storage  space 
is  desired. 

K(12) , L(12) : 1(35)  and  1(37)  are  constants,  A1  + 1 and 

1,  respectively,  used  as  stepping  indices  for  the  common 
region  of  (S) . 1(35)  is  the  interval  between  elements  in 

successive  rows  of  the  same  column,  and  1(37)  between 
elements  in  successive  columns  of  the  same  row.  Note  that 
1(36),  on  the  spine  of  (I),  is  the  locator  for  the  common 
region,  here  set  arbitrarily  at  831. 

K(13),  L (13)  : 1(39)  and  1(40)  are  constants,  A1  + 1 and 

Dl  + 1,  respectively.  1(38)  is  the  interval  between  elements 
in  successive  rows  of  the  same  column  of  the  a region  of  (S) 
and  1(40)  is  the  interval  between  elements  of  successive 
columns  of  the  same  row  of  the  6 region.  These  two  are 
interchanged  when  (S)  is  transposed.  They  are  used  as 
stepping  indices. 

K ( 1 4 ) , L(14)  , K ( 15) , L (1 5 ) : 1 (41),  1 (43),  1 (44),  and  1(46) 

are  not  used  regularly  in  this  version  of  the  program.  How- 
ever, the  spaces  1(41)  and  1(43)  are  used  to  store  strategy 
locators  needed  in  SUB  4000.  The  spine  of  (I)  is  empty 
except  for  1(27)  and  1(36),  which  have  already  been  discussed. 


45 


and  Tor  1(1)  and  1(3)  which  are  used  as  working  storage  for 
step  indices  in  the  pivot  region  of  S.  These  vary  from 
target  group  to  target  group.  In  SUB  8760,  one  of  the  pair 
is  set  at  1 and  the  other  at  N(G9  + G)  , these  being  the 
column -to-coluron  and  row-to-row  steps  for  target  group  G, 
as  shown  in  Figure  8. 
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This  loop  reads  data  on  the  number  of  targets  T(G)  and  the 
value  of  each  target  V(G).  It  initializes  Q(G)  at  '.ne  value 
V (G) . The  loop  computes  VI  = IV(G)  as  the  simplest  way  of 
making  VI  > max  (V(G) } for  use  in  the  testing  processes. 
Finally,  it  initializes  the  < G)  locators  to  the  (S)  matrix, 
whose  current  dimensions  are  zero  because  there  are  no  excess 
elementary  strategies  and  no  weapons  to  start  with. 
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This  doable  loop  reads  the  matrix  of  kill  probabilities  into 
(P)  and  prints  it  for  reference. 
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This  loop  reads  the  intercept  probabilities,  prints  them 
for  reference,  and  stores  the  complements  in  (D) . 

C.  Control  Section 

This  control  routine  directs  the  computations  over  a 
path  segment  and  is  repeated  for  each  segment  of  the  path. 
It  consists  of  four  parts: 

Segment  initialization  (1000-1380) 

Strategy  on  boundary  (1500-1590) 

Strategy  in  region  (1700-1810) 

Segment  termination  (1820-1880) . 

The  two  middle  parts  are  repeated  in  alternation  as  many 
times  as  required  along  the  segment. 
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The  path  segment  is  defined  by  two  numbers  in  this  version 
of  the  program:  PI,  the  identifying  number  of  the  single 

weapon  type  to  be  varied  as  parameter  along  the  segment. 
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and  P3,  the  terminal  number  of  weapons  of  that  type.  Here, 
the  numbers  are  input  by  the  operator,  but  the  program  can 
be  modified  to  read  the  numbers  as  data  or  get  them  from  a 
data  file.  Meaningful  values  of  PI  are  integers  from  -Di  to 
+A1,  with  defense  weapon  types  identified  by  a - sign  f-nd 
the  end  of  the  program  run  signaled  by  a zero.  P3  may  have 
any  non-negative  value.  Thus,  the  number  of  type  PI  weapons 
may  be  increased  or  decreased  along  a path  segment.  However, 
a decrease  to  0 should  be  avoided  because  of  the  possibility 
that  round-off  errors  will  throw  the  path  into  negative 
regions  of  the  resource  space.  Special  protections  against 
this  accident  were  written  into  an  earlier  version  of  the 
program,  but  have  been  eliminated  in  this  version.  If 
Pi  = 0,  the  run  stops;  otherwise  it  continues. 

U4f  s 1 = M ) 

I Cb  C 

lCt(.  (.DSU* 

1C  VC  3- 5'C  *<:<>♦♦  1)  ) 

1C7S  l-h?) 

if h r Usapsc P3-5<k< > > 

\C9(.  ih  t*<»r  HfctN  ii  ri o 

noc  bi-iUh  *pvr 

HIT  WiTC  HK 

This  sequence  sets  Si  at  +1  or  -1,  converts  Pi  to  a positive 
number,  and  goes  to  SUB  8690  to  set  the  (K)  and  (L)  indices 
associated  with  SI.  It  also  sets  P2  at  +1,  0,  or  -1  depending 
on  whether  the  number  of  weapons  is  to  be  increased,  unchanged, 
or  decreased  from  the  number  at  the  end  of  the  preceding 
path  segment;  sets  V2  at  0,  VI,  or  2V1  for  later  use  in 
testing;  and  sets  P4  at  the  amount  of  change  in  the  number 
of  weapons.  If  P4  = 0,  the  program  goes  to  SUB  9290  to 
compute  a strategy,  and  then  branches  to  the  segment  termina- 
tion, an  option  used  when  the  operator  has  just  completed  a 
path  segment,  obtaining  a printout  of  either  the  attack  or 


defense  solution  strategy,  and  wishes  to  obtain  a printout 
of  the  other  solution  strategy  at  the  same  point  of  the 
path.  Otherwise,  the  run  continues. 
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This  sequence  examines  the  proper  weapon-number  row  of  (A) 
tc  f\.nd  the  proper  column  for  weapon  PI,  designating  this 
column  by  R5.  If  weapon  PI  already  has  a column  specified, 
initialization  is  completed,  and  the  program  branches  to 
1710.  If  not,  cher.  the  column  is  picked  that  will  put  PI 
in  the  proper  sequence  with  the  weapons  already  listed  in 
(A).  For  example,  if  attack  weapon  1 is  to  be  added  to 
Figure  6,  R5  will  be  1,  and  the  program  will  continue  at 
1180. 
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This  sequence  opens  up  and  zeros  the  R5  column  of  (S) , 
shifting  other  columns  as  required.  Statements  1200-1250 
open  a column  (or  row)  in  the  common  region,  and  statements 
1254-1264  open  a column  (or  row)  in  the  ot  (or  6)  region, 
depending  on  whether  Si  = +1  or  -1. 
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This  sequence  opens  up  and  zeros  the  R5  column  of  (A) , and 

enters  the  number  PI  in  the  weapon  row  of  this  column. 
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This  sequence  steps  two  indices  in  (1)  to  reflect  the 
addition  of  a weapon  type,  and  also  steps  the  variable  Rl, 
which  measures  the  greater  of  the  two  quantities: 

1(2)  + Excess  Attack  Strategies 

1 (4)  + Excess  Defense  Strategies. 

Rl  is  automatically  initialized  at  zero  and  then  is  increased 
or  decreased  as  required.  Rl  is  a measure  of  the  size  of 
(S) , and  is  used  in  setting  the  current  size  of  the  third- 
stage  matrix  (U) . 
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This  sequence  controls  the  computation  of  strategy  crossing 
a boundary  (see  Section  III.D.),  where  the  Q-basis  is 
rectangular  and  the  parameter  is  marginal  value,  as  indicated 
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by  the  flag  S2  - -1.  The  flag  SI  ® 1 Indicates  an  attack 
strategy,  SI  - -1  a defense  strategy.  SUB  9000  computes  (X) 
and  (¥) . SUB  2000  finds  two  candidates,  H2  and  H3 , for  a 
critical  value  of  the  parameter.  These  are  compared  and 
the  program  branches  to  SUB  6000  to  add  a row  (or  column)  to 
the  basis,  or  to  SUB  7000  to  delete  a column  (or  row)  from 
the  basis.  In  the  latter  case,  R1  is  reduced  because  (S) 
has  changed  size.  In  either  case,  the  new  basis  is  square 
and  the  program  continues  at  line  1700. 
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This  sequence  controls  the  computation  of  a strategy  in  the 
in  terior  of  a region  or  at  a regional  boundary  (see  Sec  cion 
III.C.),  where  the  Q-basis  is  square  and  the  parameter  is 
resource  variation,  as  indicated  by  tne  flag  S2  = 1.  The 
control  sequence  is  the  same  as  in  the  1500  routine.  However, 
H2  and  H3  are  initialized  at  P4  instead  of  at  V2,  R\  is 
increased  if  a row  (or  column)  is  added  to  the  basis,  and 
a test  at  line  1790  provides  an  exit  to  the  sequence  for 
terminating  the  path  segment.  Otherwise,  the  program  con- 
tinues at  line  1500,  which  alternates  with  1700  until  the 
pa tli  segment  is  terminated. 
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The  path  segment  is  terminated  by  posting  the  current  force 
level  in  ($),  adjusting  a right-side  element  of  (S) , and 
printing  the  terminal  strategy  in  SUB  5000.  The  program 
then  goes  back  to  line  1000  for  the  next  path  input. 

D.  Test  Subroutine 

The  test  subroutine  finds  two  candidates  for  the  critical 
value  of  a parameter,  as  discussed  in  Section  III.  Candidate 
values,  H2  and  H3,  are  initialized  in  the  control  program 
and  modified  during  testing.  H2  is  associated  with  the  non- 
negativity condition  and  H3  with  the  scalar-product  condition. 
The  final  choice  between  candidates  is  made  after  return  to 
the  control  routine. 

The  subroutine  consists  of  <_n  index  section,  a branch 
contro1 ling  tests  on  a defense  strategy,  and  a branch  con- 
trolling tests  on  an  attack  strategy.  Seven  auxiliary 
subroutines  serve  both  control  branches. 
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This  sequence  initializes  indices  for  the  first  target 
group:  II  locates  the  intercept  elements  in  (X)  and  (Y) ; 

JO  locates  the  control  row  of  (A) ; J1  locates  the  base 
elementary  strategy  in  {A) ; and  14  locates  the  row  of  (P) 
containing  single  shot  probabilities  of  kill  for  targets 
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of  this  group.  In  addition#  1(42  + Si)  is  set  at  0 , and  J7 
is  identified  with  42  - SI  so  that  I(J7)  may  be  used  as  a 
working  strategy  locator  in  connection  with  the  value 
function,  SUB  4000.  After  initialization,  the  program 
branches  to  one  or  the  other  of  the  teot-control  sequences. 
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This  sequence  controls  the  testing  when  (X)  and  (Y)  repre- 
sent a defense  strategy.  The  test  sequence  for  each  target 
group  is  illustrated  by  the  flow  diagram  of  Figure  11.  SUB 
3000  conducts  the  non-negativity  tests  for  components  of 
the  defense  strategy.  The  rest  of  the  diagram  is  concerned 
with  the  scalar-product  tests  for  elementary  attack  strategies 
not  in  the  basis. 

The  general  scheme  is  to  generate  an  elementary  attack 
strategy,  compute  the  value  function  for  its  combination 
with  each  elementary  defense  strategy  in  the  basis,  compute 
the  scalar  product  with  (X)  and  (Y) , and  test  for  criticality. 
The  generating  scheme  in  this  version  of  the  program  is  what 
we  call  a "floating  lid."  It  generates  all  strategies  over 
a truncated  rectangular  region  whose  size  is  controlled  by 
a test  control  row  of  the  (A)  matrix.  In  the  example  of 
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Test  Criticality 
SUB  3300 


FIGURE  11 

TEST  SEQUENCE  (Statement  Lines  2110-2240) 
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Figure  7,  the  test  control  row  for  G = 3 consists  of  the 
numbers  2,  2,  2,  4.  All  test  strategies  are  generated  that 
would  not  increase  any  o £ those  numbers  by  more  than  1, 
specifically  all  strategies  satisfying: 

0 : a,  < 3 , 

0 < a,  < 3 , 

0 < a5  < J , 

0 < la  < 5 . 

“ m 
m 

The  first:  three  conditions  define  a rectangular  region  of  64 
strategies,  and  the  fourth  truncates  it  so  that  only  44  test 
strategies  are  actually  generated.  The  "floating  lid1’  scheme 
is  a compromise  that  gives  excellent  results  for  the  two-sided 
game  but  may  be  slightly  off  optimum  for  one-sided  attack 
allocations,  especially  if  some  of  the  weapon  types  have  very 
low  kill  probabilities. 

Test  attack  strategies  are  generated  recursively  in  the 

1-30  region  of  matrix  (A) , and  the  current  critical  candidate 

is  stored  in  the  31-60  region  of  (A).  The  order  of  storage 

is  a values,  la  , and  computed  v's,  one  v for  each  elemen- 
m m m 

tary  def-nse  strategy  on  the  target  group  under  test.  The 
method  of  generation  is  by  a nest  of  implicit  loops  on  the 
weapon  types  currently  in  play. 

The  scalar  product  of  the  test-strategy  vector  with  the 
solution  vectors  is  represented  as  X2  + H*Y2,  where  the 
coefficients,  X2  and  Y2,  correspond  with  (X)  and  (Y) . For 
efficiency,  partial  sums  are  computed  by  recursion  as  XI 
and  Yl. 
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We  can  now  describe  the  flow  diagram  of  Figure  11  in  more 
detail.  SUB  3200  generates  an  initial  test  strategy  with  all 
zero  elements.  SUB  3700  computes  values.  If  Y2  i,  0,  the 
strategy  cannot  be  critical;  otherwise,  SUB  3300  compares  it 
for  criticality.  The  A(K(10))  test  determines  if  the  "lid" 
on  the  sum  has  been  reached.  If  not,  the  first  weapon  type 
is  set  by  W * 1.  If  the  lid  has  been  reached  and  W < K(l), 
then  SUB  3500  reduces  A(W)  to  zero  and  sets  W = W + 1.  The 
A(W)  test  determines  if  the  "lid"  on  weapon  W has  been 
reached.  If  not,  SUB  3600  steps  A(W)  and  the  cycle  repeats. 

If  so,  and  if  W = K(l),  then  testing  on  the  group  is  finished, 
and  SUB  3820  steps  indexes  for  the  next  group. 
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This  sequence  controls  the  testing  process  when  an  attack 
strategy  is  being  determined.  It  is  evident  by  comparison 
with  (2100-2260)  that  the  general  scheme  and  computational 
subroutines  are  the  same,  but  there  are  significant  differ- 
ences in  the  control  processes.  If  no  defense  weapon  has 
been  introduced,  if  G is  an  undefended  group,  or  if  G has 
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not  been  brought  under  attack*  no  tests  of  new  defense 
strategies  are  appropriate.  So*  an  immediate  branch  to 
2660  occurs. 

The  range  of  elementary  defense  strategies  to  be  tested 
is  determined  solely  by  the  attack-summation  test  control 
number  A (JO  + 1(31)),  and  the  sole  cut-off  control  is  state- 
ment 2640  which  prevents  defense  use  of  more  weapons  than 
A (JO  + 1(31)).  In  the  example  of  Figure  7,  defense  strate- 
gies to  be  tested  must  satisfy  the  single  condition: 

0 < £6  < 4 . 

n n 

If  there  were  two  types  of  defense  weapons*  this  inequality 
would  call  for  the  generation  of  15  elementary  strategies. 

E.  Auxiliary  Subroutines 

3U(  K/,s\(  h ( 7 ) ♦ u> 

3 f 1 ( fF  i H 3 1 f f 

3<  I = I 1*1  V.  1 1 + « 4 + ) 

3 ( 3 ( If-  f ( I ) > - . c ; ( l Ht.'J  3 C **  ( 

3f  4C  H--XCI  >/r'C  I ) 

3 ( b 0 I I-  H r < H 
3 ( e l 3 ;•  = H 
3 < 'i<  ( psi 

3(  '-if  f - 1-  I l 
3(*r  \im  I r 
■*  1 C C ^ fc  1 u f \ 

This  subroutine  tests  the  condition  x?  2 0 for  elementary- 
strategies  in  the  basis.  It  solves  for  H,  the  value  of  the 
parameter  that  makes  x?  = 0*  and  records  H2,  the  candidate 
for  critical  value.  It  also  records  G2,  the  group  in  which 
that  candidate  is  found,  and  M2,  a locator  for  the  particular 
elementary  strategy  among  those  of  group  G2. 
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ns* ft*  v. - 1 in  k f if  > 

3MC  Affc)=C 
3??C  VfcXI  V. 

?y*C  K\tK'rS/(  I I ) 

3ibC-  r 1,  T'J-sYf  I I ) 
nj'fiC  Mr  I UhlNi 

This  subroutine  generates  the  null  strategy  as  an  initial 
test  strategy  and  initializes  the  scalar  products. 


330( 

an ic  t h *<3<h  ihfn  3ao( 

3 3P(  H3  = H 

aaai1  i’3st,i 

337f  K?m  K s 1 lk  K{|f)*t<A»| 

33bC  A<K+3C)s^<K) 

339 f \lfcXT  K 
34(-(  hfc.iU/v>J 

This  subroutine  computes  the  value  of  the  parameter  that 
makes  the  scalar  product  zero,  compares  this  value  with  H3 
and  replaces  H3  if  H3  i H;  irt  this  case,  the  subroutine 
records  G3  - G and  transfers  all  pertinent  data  from  the 
1-30  region  of  (A)  to  the  31-60  region,  where  the  best  test 
strategy  is  recorded.  This  completes  the  testing  of  one 
elementary  strategy. 


3srr  x i>x  i-Af  v.)*xor.) 

3 b l c r i = r i-a< v »*y< y > 

3b?0  A<K<  1C>  >*ACK<  ID  )-A<fc> 

3b3r  A(V.)=t 
3bz,c  V = V + 1 
3 b b C utiUhiN 

This  subroutine  reduces  A(W)  to  zero  after  making  related 
changes  and  steps  W to  W + 1. 
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at  tc  <-(k<  ir»>st(K<  ic>  > • i 

C A Ii^s<|«(  ..) 
y i»  y f=y  i ♦ y*  v ) 

3 r *. C i'(  ium 

This  subroutine  generates  a new  test  strategy  by  stepping 
A ( W)  and  making  related  changes. 

;*/(»  «*|| 

:»/«(  i (//)-.  1 1 

,'i  v ,'u  y i . n x - i i . . ■<*)*  l 

a * •.  t y * i 

j , \ t i c .1  / )=  i < J / 5 *“(  v ) 

,V/ <•  C WOi 

awe  <■(«(  io*k>:v 

y i ■<  r \ ; - = < > ♦ v * X ( i*  ) 

• / •* . r = r •,  ♦ \,  » y ( * ) 
c i a /.  i 
a < i • it-  i u n.s 

This  subroutine  controls  access  to  SUB  4000,  the  value 
function;  it  also  sums  for  the  coefficients,  X2  and  Y2, 
of  the  scalar  product. 

'!  - r ( i I = I I ♦ •<  *.  ♦ ; • 

t .li  =,):.*  I c ) 

. t i = . J I ♦ * ( V ) * ' •<  A * i ) 

:wtc  i a- i 4«v i 

f fklu.-N 

This  subroutine  steps  indices  for  the  next  target  group. 

F.  Value  Subroutine 

This  subroutine  is  addressed  from  statement  3760  only. 
It  computes  and  returns  one  number,  V,  the  expectel  value 
survivinj  of  a single  target  of  type  G when  it  is  attacked 
using  the  elementary  strategy  located  at  1(41)  in  (A)  and 
defended  using  the  elementary  strategy  located  at  1(43)  in 
(A) , these  two  locations  having  been  defined  at  statements 
2C30  and  3750,  respectively. 
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In  this  version  of  the  program,  it  is  assumed  that  the 
firing  sequence  of  attack  weapons  is  worst  to  best  with  the 
idea  of  using  poor  weapons  to  exhaust  interceptors , and  the 
sequence  of  defense  weapons  is  best  to  worst  with  the  idea 
that  any  leftover  interceptors  will  be  the  worst.  This  is 
the  order  in  which  weapons  have  been  numbered  in  the  original 
data  input  for  the  program.  Interceptors  and  attacking 
weapons  are  thus  matched  in  pairs,  as  is  illustrated  in 
Figure  12  for  the  strategies: 

a = (3,2,3) 
t = (2,4)  . 

In  the  illustration,  the  expected  value  surviving  is: 

V * V 
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where  PKT^  is  the  probability  of  target  kill  by  a single 
attack  weapon  of  type  i and  PI^  is  the  probability  of  inter- 
cept by  a single  defense  weapon  of  type  j . Subroutine  4000 
computes  V according  to  this  scheme. 


A ( 'C  (,  V = V ( b ) 

«c ir  o, i =( 

Ai?C:  K'H  P=  J IC  IC4) 

A C 3 C A7  = £M(4l)+*0 
ACIO  It-  £7  = ( l-tfrN  A's 7C 
I b= I A+£(  I C ?3> ♦ A) 
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FIGURE  12 

A PAIR  OF  ELEMENTARY  STRATEGIES 


This  sequence  initializes  value,  defense  weapon  type,  and 
number  of  interceptors;  begins  a loop  on  A extending  to  line 
4270;  set3  A?  = the  number  of  type  A weapons  in  the  a 
strategy;  if  A7  = 0,  goes  to  the  next  A;  otherwise,  sets 
15  to  locate  PKT  for  weapon  type  A in  matrix  (P) . 


4C-9C  I*  OP.*  T8i.v  4Pt( 

4 ICC  IP  l>C  7 8£.\  4 1 fit 
4 1 1 C IP  1'slCSO  7HfcN  4?tC 
4 1?C  U=L>*  1 
4130  I a A(  I < 43 ) ♦ L') 

4 1 V C t»U  T fc,  4 ICC 

This  sequence  decides  whether  an  interceptor  reading  is  in 
order;  if  the  target  group  is  undefended,  the  sequence 
branches  to  4260  to  compute;  if  I > 0,  it  branches  to  4180 
to  test  A7;  if  the  defense  is  exhausted,  it  branches  to 
4260  to  compute.  The  sequence  also  steps  D,  sets  I = the 

-fr- 

number  of  type  D interceptors  in  the  6 strategy,  and  returns 
to  4100  to  see  if  I > 0. 


4 18  C It-  A 7 < = I I8fc.NI  4f-lC 

4 19  0 J=I 

4fcOC  GCTO  4?rjC 

4?1C  J=A7 

4PPC  A 7 = A7 - J 

4P3C  I=l-J 

4 ? 4 C t=V/4(l-p(Ib)*L,PA(l(?b)+ti)))»J 
4^bC  COTS;  4C7C 

This  sequence  sets  J = min{A7,I};  reduces  both  by  J;  computes 
a new  value  of  V at  4240;  and  goes  back  to  4070  to  see  if 
A7  = 0. 

4 ? 6 C-  9=W(  1-PCIb)  )»  A7 
4 ?7  C \LX I A 
4^8  C M-.  JO/mM 
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This  sequence  computes  a new  value  of  V when  there  is  no 
more  defense.  When  all  of  the  attack  weapons,  but  not 
necessarily  all  of  the  interceptors,  have  been  exhausted, 
the  program  returns  the  computed  V to  statement  3770. 

Some  general  observations  may  be  useful.  SUB  4000 
contains  the  only  two  statements  in  the  entire  program 
that  make  direct  use  of  the  probability  data  stored  in 
matrices  (P)  and  (D) . This  feature  gives  the  operator 
some  latitude  to  simplify  4240  and  4260  by  precomputing 
some  of  the  factors  and  storing  them  in  (P)  and/or  (D)  at 
the  time  of  initialization.  Thus,  he  might  save  computer 
processing  time  by  using  more  storage  space  and  some  extra 
indexing.  In  fact,  the  current  version  of  the  program 
contains  at  statement  540  a precomputation  of  the  complement 
of  PI,  but  that  doesn't  require  any  extra  storage.  We  felt 
it  more  desirable  at  this  time  to  minimize  multidimensional 
storage  requirements  in  (P)  at  the  expense  of  added  pro- 
cessing time.  However,  if  a great  many  runs  were  to  be  made 
on  small-scale  problems,  it  might  be  judged  worthwhile  to 
increase  the  amount  of  precomputation  in  the  value  function. 

Of  much  greater  importance,  however,  is  the  flexibility 
that  the  operator  has  to  use  an  entirely  different  value 
function,  even  going  sc,  far  as  to  read  in  all  the  values  as 
arbitrary  input  data.  In  case  a different  value  function  is 
to  be  used,  the  operator  should  provide  the  following  general 
modifications  in  the  program; 

-1  Change  the  initialization  procedures  to  read  in 

the  desired  input  data  and  to  precompute  and  store 
the  desired  quantities  in  (P) , (D) , and  any  other 
unused  matrix. 
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-2  Replace  SUB  4000  with  a subroutine  designed  to  get 
V for  any  pair  of  strategies  located  by  I (41)  and 
1(43)  at  statement  3760.  Be  careful  to  avoid 
accidently  changing  any  of  the  variables  used 
elsewhere  in  the  program.  In  general,  the 
undifferentiated  letter  variables  are  available 
for  use,  with  the  exception  of  G,  W,  K,  and  M, 
which  have  specified  values  at  the  time  of  access 
to  4000.  Obviously,  A7  and  IS  are  also  available. 
Other  variables  should  be  used  only  after  care- 
fully examining  the  entile  program  for  conflicts. 

-3  Make  whatever  indexing  changes  are  consistent 

with  the  new  value  function;  specifically  change 
the  setting  of  14  at  2030  and  its  stepping  at  3850, 
if  desired. 

-4  No  other  changes  need  he  made  as  long  as  the 
operator  is  satisfied  with  the  "floating  lid" 
method  of  generating  test  strategies.  That  can 
be  changed  too,  but  this  is  not  the  place  to 
discuss  how. 

G.  Print  Subroutine 

In  this  version  of  the  program,  the  print  subroutine  is 
addressed  only  from  statement  1870  at  the  end  of  a segment. 

It  prints  either  an  attack  strategy  or  a defense  strategy, 
as  determined  by  the  current  value  of  SI. 

5CCC  L'EF  FNA(  J )-H*Y<  I >*X<  1 ) 

This  line  defines  FNA(I)  as  a function  to  be  used  in  combining 
the  (X)  and  (Y)  solution  vectors  into  a single  vector  by 
means  of  H,  the  value  of  the  parameter. 
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b<  V = (I 

This  line  initializes  the  value  of  the  game. 

bCPC  Uh  U1  it  K(l)  _ 

bC3C  r~b 

bC*C  HsA(L(K)>k) 

SCbC 

b{;<LC  J-JvlN  1 bi>). 

b C.  V 0 i\  k X 1 l. 

This  sequence  computes  the  marginal  value  for  each  of  the 
opposition  weapon  types  currently  in  play;  adds  the  product 
of  marginal  value  by  number  of  weapons  to  V;  and  prints 
weapon  identifying  number  and  marginal  value. 


t ch c i =k(  n 

5 09  C ,J  = K(H)+K(0)*(,| 

b I r f h <fj  I'  l>-  1 1 \j  I;  | 

This  sequence  begins  a loop,  ending  at  statement  5210,  that 
computes  and  prints  the  augmented  strategy  for  each  target 
group. 


b 11 C 1=1+1 

b 1 P t M(Otl'i'iMl) 

b 1 3C 

b 1 4 C P /->  I .vi  1 J l'»  *■/ 

This  sequence  computes  the  negative  of  the  per  target  inter- 
cept value,  multiplies  by  T(G)  to  get  the  total  intercept 
value  for  the  group,  adds  this  to  V,  and  prints  the  group 
number  and  the  group  intercept  value. 


5150  Hih  N*0  T0  N < K ( 7 > ♦ G) 
bltC  I « I ♦ 1 
5 170  J* J*K  < 9 > 

517b  f'KlNT  li 
5 180  H4h  L * 1 TtJ  Lt  1 ) 

5185  f-MNl  M J + L > J 
5190  NEAT  L 

5195  PKINT  T<G>*ENA< 1 > 

5 POO  NEXT  N 
5 P 1 0 NEXT  b 


For  each  elementary  strategy  on  a target  group#  this  sequence 
prints  the  number  of  each  type  of  weapon  per  target#  and  the 
number  of  targets  on  which  that  strategy  is  used. 


5PPC  HKINT  “V5=MV 
5P30  HETUWVJ 

Line  5220  prints  VS#  the  value  of  the  game. 

Section  VI. B.  contains  many  illustrative  printouts. 

H . Addition  Subroutine 

This  subroutine  is  addressed  from  statement  1560  or 
statement  1770  when  the  critical  value  of  the  parameter  is 
H3  (implying  that  a new  elementary  strategy  must  be  added 
to  the  Q-basis).  It  makes  the  necessary  changes  in  (A),  (R) , 
(S) , and  the  associated  indices. 

The  target  group  affected  is  G3,  as  recorded  at  3330, 
and  the  numerical  elements  of  the  new  strategy  are  on  hand 
in  the  31-60  region  of  (A) . 


<000  G=  63 
fcCJC  S 1 = - £>  J 
fcOPO  G0SUB  869C 
6C3C  ODSUL'  R76C 
6C3S  GCbUE*  H835 
60*0  Sls-Sl 

This  sequence  identifies  the  target  group  and  gets  the  correct 
indices  from  subroutines  8690,  8760,  and  8835.  For  this 
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purpose,  SI  is  temporarily  reversed  to  orient  (K)  and  (L) 
properly  for  adding  either  a row  or  a column.  The  reader 
will  find  it  helpful  to  visualize  the  subroutine  as  adding 
a row,  and  the  text  will  follow  this  line  of  thought,  with 
a few  identified  exceptions. 


<>rsc  j?=j  i*Kt *? > * c K4*  n 

#CfO  Kft  JsJV  lu  JJ-'+l  STfcf-1 

iCi a a<j*k<9>>  = m.j> 

t ( H (.  ,\l  fc  X 1 J 

H.vG  f U K = I It'  L ( in 

f icr  pc  jp+n  > = 3 on  > 

6lir  I h A(  JON)>f»<  M *\ > 1WFN  f!3C 
MPC  JON  > = Af  ION  ) 
f.  1 3 ( N fc  X 1 N 

This  sequence  makes  all  necessary  changes  in  (A)  and  locates 
the  base  strategy  row  of  the  G3  + 1 target  group  in  (A) . 

This  is  the  row  where  the  new  elementary  strategy  being 
added  to  the  G3  group  will  be  stored.  The  program  makes 
room  for  the  new  strategy  by  shifting  upward  all  higher  rows. 
It  stores  the  new  elementary  strategy  and  the  weapon  sum  in 
the  J2  row  and  changes  the  test  control  elements  in  the  JO 
row  as  required  by  the  "floating  lid"  method  of  controlling 
tests.  (A  different  method  of  controlling  tests  might  call 
for  modifying  these  operations  on  the  JO  row,  or  even 
removing  them  from  the  program. ) 


<lbC  hluh  «a(  <K(  (.)  *( 9 ) - j If)  (-(  K c l ) MFf-  -I 

f 1 b t.  ht  tf  + K ( I3))si-(M) 

(.  Iff)  oO  + K(  13)  ) = SCI*) 
t 17C  Ntxi  y. 

This  sequence  shifts  elements  of  the  a region  of  (R)  and  (S) 
to  open  the  proper  row  for  storing  the  new  strategy. 
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fclrtC  hdh  I«6C(.9»-t  (»<(♦!)  bTEt-l 

t |Hb  M 1*L4)»M  1 ) 

6 19 C S< I*L4)=S< 1 ) 

6 |9  S (MtXI  1 

6200  If  :>1  = -1  THfc \ t.290 

This  sequence  shifts  elements  in  the  pivot  regions  of  (R) 
and  (S)  to  open  a row  for  storing  the  new  strategy,  performing 
the  shift  for  those  pivot  regions  pertaining  to  groups  with 
G > G3.  If  a true  row  is  to  be  added,  as  indicated  by 
SI  * -1,  the  program  branches  to  6290  since  the  required 
space  is  now  open.  This  case  can  be  illustrated  by  Figure  8. 
When  elements  7-12  are  shifted  to  10-15,  spaces  7,  8,  and  9 
become  available  for  a new  row  of  the  first  group. 


*2Cb  K3)*I<3>  + 1 

6Pic  i=c;cb+»> 

6J1S  J=L* 

622C  hiirt  L = ? Iv,  L 4 
f.f-22  J = J-  1 

Kh  K=|  lb  K 4 
f 23C  1=1-1 
t?4Z,  h(  I*J>=H<  I ) 
t ?bC  1 ♦ J ) =i< 1 ) 

6 2 6 C NEXT  X 
b?HC  N EX  1 L 

If  a true  column  is  to  be  added,  as  indicated  by  SI  s 1,  this 
sequence  performs  variable  shifts  to  open  the  proper  spaces 
within  the  G3  pivot  region.  As  an  example,  consider  the 
addition  of  a column  to  the  first  group.  The  elements  7-12 
are  shifted  to  9-14  by  the  6180  sequence.  Then  the  6205 
sequence  shifts  4-6  to  5-7,  so  that  spaces  4 and  8 become 
available  for  the  new  column.  The  index  1(3)  is  stepped  to 
reflect  the  increased  row-to-row  interval. 


629  (.  *SMt(K  <£)♦(>♦  I) 

6 3CC  R(rib)*C.-{  i.3)-£C31«l<  1(<>  > 
fc  3 1 C MMb)  = T< b3>*MMS> 

633C  hk/h  L:  1 Tk!  LC  1 ) 

63^0  M Mb»L  ) = T<  u3>*<  #>(  3(HL  ) -£(  J 1*L  ) ) 

63bC  NfcX7  L 

This  sequence  enters  the  first  (or  right-side)  element  of  the 
newly  opened  row  of  the  a region  of  (R)  in  the  form  vt 1 - v^ 
without  multiplying  by  T.  It  then  enters  the  same  element 
of  (S) . In  this  version  of  the  program,  all  elements  of  (S) 
contain  the  scaling  factor  T(G)  , but  none  of  the  (Q)  or  (R) 
elements  are  scaled.  This  sequence  enters  the  scaled  a 
differences  in  (S)  . These  differences  are  not  saved  in  (R). 


t3fcC  ft*  (.=  G3+  ] !«;  (yi 

6 36  b tf<  k ( 7 ) ♦(•)  = «<  K ( 7 ) ♦ C-)  ♦ I 

63  7 C GCK<  r ) *l)*K<  J 3 ) 

63  Vb  0(  l»J  =CC  0)  *L4 
6 3 H C.  Nk/1  i, 

6 39  C N ( K ( 7 ) + ( 3 > =N  < K ( 7 ) ♦(-:  3 ) ♦ J 
K.C.C  If-  C3>U4  HF,\  67 3 C 

This  sequence  readjusts  all  indices  affected  by  the  addition 
of  a row  to  G3.  If  G3  > G4,  the  program  branches  to  a RETURN 
statement,  since  no  further  computations  are  needed  on  an 
undefended  target  group,  which  has  no  excess  defense  strategies 
and,  hence,  no  pivot  region  in  (S) . 


64f.b  K4=K4*  i 

<4  K K-K(b)+iv!(K(7)+(>9)  li.  Kf+K  4 ilFt  -) 

( 4P0  i‘.  CK  + l)=.iCK  ) 

64 bf-  iNFXT  K 
t AtC  aCKOK4>  = f 
647  r.  Ms\c 

648C  I,  I b=I 0+ I ( 1 > *<K4-  1 ) 

649C  J=  3 1 +L ( »C) 


65CC  b0h  L * 1 TU  LA 
6 53C  .Is,/*  I 

6b AC.  hC  I >sA<J>-£<31*L<  1C)  >*KCNJ  ) 
t 5*5  3C I > s7  Cb3)*K<  » > 

6550  1 « I ♦ I C 3> 

£555  M|s,n|*L(  1 J) 

65fcC  NEAT  L 

This  sequence  sets  K4  = the  new  number  of  rows  in  the  pivot 
region  and  readjusts  the  pivot  index  (0).  It  enters  new 
values  in  the  open  row  of  the  pivot  regions  of  (R)  and  (S) . 


fcb/O  L6=G 
£5  75  L'/s  l 
65RC  L3=l<3) 

6585  I £ = I 5 

6590  bAb  L I = 1 Tt:  La 

1 6 1 c n o<lc*ld  = o then  <r-u 

&t?C  hJ  = H(LC*L1) 

l 63C  L = T((i3)*K(  l',‘) 

ttAC  3(1 £)  = ;><  IP) -6 

IbAb  M3  = MC*K(  )3>*(X 1-  J > 

t 6 5C  I 3= I 0* 1 C 1 ) *(K  l-  1 ) 

( tfcC  UoiUl-  RS3P 
667  0 l.UJtt  6 6S<r, 

< 68 C L£=L  1 
EfcVC  I pa  I P* I ( 3) 

6695  NEXT  LI 

This  sequence  operates  on  the  new  row  of  (S)  to  reflect  the 
effect  of  pre-existing  pivots  in  the  other  rows  of  (S).  It 
initializes  L6,  a variable  used  to  record  any  column  without 
a pivot  that  may  be  found  during  the  ensuing  process;  pre- 
sets L7  and  L8,  interval  step  indicators  needed  in  SUB  8532; 
and  initializes  12  as  the  first  element  of  the  new  row  in 
the  pivot  region.  It  loops  on  columns.  If  there  is  no 
pivot  in  the  LI  column,  the  program  branches  to  6680, 
records  L6  = LI,  steps  12,  and  goes  to  the  next  column.  If 
there  is  a pivot  in  the  LI  column,  the  program  records  its 
row  as  Kl,  presets  Q,  a multiplier  used  in  SUB  8532,  replaces 
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S (12)  by  S(I2)  - Q,  locates  the  first  elements  of  the  K1  row 
in  the  a region  and  the  pivot  region , goes  to  SUB  6532,  and 
then  goes  to  the  next  column. 

t 7CC  I ¥ L 6=  C iHfc.N  t /3t 
<7C5  K s k 4. 

<*• v i r>  i = 

» t'c  C If  AhSCtC  I > >>•  ) K 1 | f 

f V 3C  nfc  J whN 

If  L6  = 0,  all  columns  have  pre-existing  pivots,  and  the 
program  branches  to  6730.  Otherwise,  it  identifies  the 
potential  pivot  row  as  K,  the  column  being  L6  > 0,  and  sets 
I as  locator  of  the  potential  pivot  element.  If  S(I)  ft  0, 
the  program  branches  to  8110  to  initiate  the  pivot;  otherwise, 
it  returns  to  the  control  section. 

I . Deletion  Subroutine 

This  subroutine  is  addressed  from  statement  1580  or  1800 
when  the  critical  value  of  the  parameter  is  H2  (implying 
that  an  elementary  strategy  must  be  deleted  from  the  hypo- 
thetical Q-basis) . It  makes  the  necessary  changes  in  (A)  , 

(Q) , (R)  , (S),  and  the  associated  indices. 

The  target  group  affected  is  G2,  and  the  elementary 
strategy  within  this  group  is  M2,  as  recorded  at  3070  and 
3080,  respectively.  The  deletion  algorithms  depend  on  the 
pivot  status  of  this  elementary  strategy. 

7CCC  b=0£ 
veil  7 t C- 

/Clb  tk/SU“  SMOb 
7 CPC  L fc  = C> 

7 C3C  JP  = J l*K(  <)}*(■;;■ 

H.4i  IF  Kr>  C.  1 M t-  iM  7bMC 

This  sequence  identifies  the  target  group  and  gets  the  correct 
indices  from  subroutines  8760  and  8835.  Since  the  (K)  and 
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(L)  indices  are  already  correctly  oriented,  there  ie  no  need 
to  go  to  SUB  8690.  The  sequence  initializes  L6,  a locator 
that  will  be  used  to  record  the  pivot  column  if  the  row 
being  deleted  is  a pivot  row  and  computes  the  location  in 
(A)  of  the  row  being  deleted.  If  M2  > 0,  the  row  being 
deleted  appears  explicitly  in  (S) , and  the  program  branches 
to  7580. 

7 CbC  I S= I C 
7 0b?  L 7 = I 
7 05  5 L«=l<3> 

7 CtC  huh  K- 1 T&  K4 
7C7C  J?*J?*K(9> 

7 09 C GOSUh  K57C 

709?  If  C:0<OK>»f  viro 

7 09  5 (j=  u-  ! 

7 ICC  Ih  CC-Cl  TMFM  71?{ 

/ I (,5  I 5= I 3+ I < 1) 

/ 1IC  \'M1  K 

If  M2  * 0,  the  base  row  must  be  replaced,  and  the  sequence 
through  line  7570  controls  the  replacement  algorithms.  The 
7050  sequence  selects  one  of  the  rows  in  (S)  to  be  the  new 
base  by  looping  on  the  rows  of  the  pivot  region,  going  to 
SUB  8570  to  compute  Q,  and  exiting  from  the  loop  when  Q ^ 0. 


7 1 ?C  Kl,.v?=K 
7130  11,13=13 

7 135  M I , M3, MbsMOK ( 13>*(H  i~  ) ) 

This  sequence  sets  locators  needed  later;  identifies  the  new 
base  row,  which  will  be  deleted  from  (S) , as  Kl,  M2;  locates 
its  first  element  in  the  pivot  region  as  II,  13;  and  locates 
its  first  element  in  the  a region  as  Ml,  M3,  M5. 
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7 ISC 

M'H  Hi  1 ll)  K4 

7 Ibb 

it  h sk  i iMhiv  7;oc 

/ J b 7 

r < M > = b ( M ) - r < k 1 > 

/ l*t 

11  = 13 

/ 1 /r 

1H  L4 

7 1 “It 

h(!>sti(!>-rClt> 

7 1 9f 

I S I ♦ I C 3 ) 

V^f'C 

1 1=1  !♦!<:*> 

O' 1C 

At.'.  1 L 

/IOC 

I 0 l s s l ^ ♦ l ( j ) 

7 JO  S 

K=M+K<  13) 

03c 

\ t At  A 

Vf.'iS 

1.  i*>e  u(  ( 7 ) - n 

7 'f 

I 1 = 13 

7 P 4 Li 

( 

7 J-M 

KiK  L * 1 I «.  (.  a 

7 > SS 

n < \ ) = l ( \ ) - f.  ( 11  ) 

7 ;- /■  ( 

1 1 = 1 1*1(3) 

7 ft  b 

\=A+I_<  13)  | 

7 J-f  / 

■Mt-A  i L 7 

This  sequence  recomputes  the  (R)  matrix  to  reflect  the  change  j 

of  base. 

and  recomputes  the  base  value  stored  in  (Q) . j 

v ; 7 c 

K.  f .*  = 1 1.'  L(  1 r ) 

7 • * ». 

A ( J i * M ) = ,-  r . K ♦ \ 1 

/ ;«<  ( 

vKl  7 

This  sequence  substitutes  the  new  base  row  in  (A) . 

/'Of 

■/  :i : 1j 

it  = i 

v c 

It  ;(-,  »*nj)  = 1'  Hr  , 7 3/1 

7 Cj  . 3 ( 

(.  ) « L f -.(■•<  * a ; ) 

v v ■ 

1 < = I 3 ♦ I ( 3 ) * ( 1 . i - i ) 

v ;< 1 . f 

Cl()=i 

/ .1/  r 

1 ; Hb.'IP 

4 

This  sequence  begins  recomputation  of  (S)  to  reflect  the 

change  of  base.  It  performs  the  computation  on  the  row  of 

(S)  identified  as  the  new  base  row.  If  this  is  a pivot  row, 

the  sequence  records  the  column  number  of  the  pivot  element 

L as  Ll , L6 

for  future  use.  Later  this  row  will  be  deleted 

f from  (S) 

but  it  must  first  be  modified  as  part  of  the 

algorithm  for  modifying  the  other  rows. 
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I'StV  iy-1  b*L<  i-  i > 

V 37b  Mbs  I ( 3t ) 

7377  L 7 sL  C 1 P) 

7373  L*sLC13> 

7 38  C FUh  K = ( T,i  W(  l> 

7 391  t*U3ui-  ib/C 
7 39  b If  K > 0 IHF.N  741C 
7 4CC  b=u-| 

/41C  .U  C(Kf*K|)  = f THf-g  7 430 

V 4<‘(.  a<  I c ) = L 
7420  OObuh  •3b3i' 

7 43b  I b= 1 b+  1 
7440  10=1^+! 

7 44b  Mb  = Mb  + K(  ) f-  > 

74bC  iNtAl  3 

This  sequence  performs  the  recomputation  for  the  rows  of 
the  common  and  6 regions. 


745b  Ibs>( 

74 Of  I ^S 1 b»I ( 2)*(L  i- 1 ) 

7 4 7*.  1-  L 7 - l 
7 4 03  L*  =1(2) 

7 4 7 b MbsV.r 

MOV  7 F 3=1  13  K 4 

7 43  0 I !•  3 = K 1 THt  .»  7 Si  0 

7 49  C U7  3U1-  H S 7 ( 

7490  IF  ; < 3 ( + 3 ) > c T-f  f m 7 s r.r 
7 49  S (.=  - 1 

7 b t C IF  ..  ( 3 C +‘*  1 ) = 0 i 4 F v 7‘.FC 
7 b u > ( !?■)  = ( 

7bW  30,'iF 

7 bl-  b I b = 1 b*  I ( 1 ) 

7 b 3 C I «‘s  1 P+  1 < l ) 

7 b 3 b M b = M b ♦ k < i ;o 
7 b 3 7 .v  F.  * 1 K 

This  sequence  performs  the  recomputation  for  the  rows  of  the 
a and  pivot  regions. 
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OA{  If  i'ltHI  *K|)  = I 1 H H <V.  /-'ll' 

/ j vj  ( . : < h t.  ♦ * i > * i ‘ ( L t ♦ f . 1 > = " 

v *.w  r i H-t-i  - \ 

/SVC  1 .1  I.:  7SL-C.- 

If  the  new  base  row  is  not  a pivot  row,  the  program  branches 
at  once  to  7820.  Otherwise,  it  sets  the  pivot  indices  to  0 
and  reduces  the  pivot  counter  P8,  since  the  row  is  going  to 
be  deleted  from  (S) . 

/•>.:  I r , C s f i H h . v / «*.  r - \ 

I r ■ l I f ♦ f : ) 

/ ' 1 . . i 7 a I l I ( i ) * ( *i{-  - 1 ) * 1 ( :*  ) * ( I-  i * I > 

,t:',  .,;is  - j 

it  Si  I !•  t :•(  I 7 ) > > . . : ! ) I He  . ■/'/%■• 

This  sequence  is  reached  only  from  7040.  If  the  row  to  be 
deleted  is  not  a pivot  row,  the  program  branches  to  7820. 
Otherwise,  it  begins  recomputing  the  (S)  matrix  prior  to 
de ' • . ng  a pivot  row,  using  an  algorithm  best  described  as 
an  '’unpivot. " The  sequence  records  the  column  number  of 
the  pivot  element  as  L6,  locates  the  pre-existing  pivot 
element,  sets  the  sensor  S3  = -1  to  indicate  an  unpivot,  and 
branches  to  7780  if  the  pivot  element  S(I9)  j-  0. 
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If  S(I9)  » 0#  this  sequence  invokes  the  unpivoting  procedure 
on  one  or  more  other'  pivots  of  the  group  in  order  to  make 
the  desired  pivot  element  3(19)  jt  0,  tests  at  7740  each 
time,  and  branches  to  7780  whenever  the  condition  is  satisfied. 
If  the  condition  cannot  be  satisfied  by  removing  all  other 
pivots  in  the  group,  there  is  an  error  printout  and  programed 
STOP  at  7770.  The  error  printout  has  never  occurred  in  any 
run  of  the  debugged  program.  However,  the  procedure  is  a 
messy  one,  and  a better  recursion  algorithm  is  needed. 


I I X ( * l - !•'  r 

/ / 1 ( I.  1 - I.  f 

i - if  In'  .-i  ■ U X 1 f f 

(••'K  , . ( X ( ♦ k | ),  i (l.  ( *1.  1 > = <■ 

If  S(I9)  jp<  0,  this  sequence  denotes  the  pivot  row  and 
column  by  K1  and  LI,  goes  to  subroutine  8160  with  S3  = -1 
for  the  unpivot  computations,  and  sets  the  pivot  index  at  0. 
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This  sequence  deletes  the  M2  row  from  (A) , (R) , and  (S) , by 
shifting  down  all  the  higher  rows.  If  SI  - -1,  the  sequence 
7090-7910  is  needed  to  close  the  separated  spaces  within 
the  G2  pivot  region  of  (R)  and  (S) . It  is  the  counterpart 
of  the  sequence  6205-6280  for  opening  spaces. 


7 7 3 b t t lu  K**  | 

v r . f k f * k > = ( k r ♦ k ♦ n 

1 | h /(Kf*k)  r(  Hl-\  V'/SS 

/ V 1 L ^(UC-M<Kl*h)>=K 
/ V t>  Ml-  A 1 >< 

7 y ( f H‘t-  Us^i»K/,  111  K(S>*v.(><rn*(9J 

/ Sty-S  |f(K)sH(l<+|) 

7 y 7 C N H X T H 

This  sequence  shifts  the  elements  of  the  pivot  index  to 
conform  to  the  new  structure  of  (S) . 


/ V-i  C K Asv<  .6-  J 

7 MO  tie  (<=(•?♦  » 7:'  I l 
7 '/  y *>  ,»•  ( * C 7 ) ♦ i ) = n ( * < 7 ) ♦ • ) - j 
« r.(. ( I ■. < K ( f. > ♦ I ) s l . ( K ( t- )♦•..)•  K ( J'<) 

-<  r c *>  i < i > = ■.  ( i • > - l ^ 

- f 1 C MM  I 0 

-i  C i f.  \ < M / ) ♦ i,f-»  s,m  < »•  ( / ) * I ; > - } 

This  sequence  changes  other  indexes  to  conform. 


i f ."it  I »•  L<  *(■  '*■  1 ( < 

H M (.  |:l(*i(3)*(l.f-|) 

^ C C I- . . r K = J 1 1 ' •<  4 

y (.7  ( 1 b .••(KC*K)>1  I ‘-1  k M 

•UoC  it-  /-KSf .sc n > » . I nt.M  iiK 
C9C  I = I ♦ I ( I ) 

Hfv1)  \ Ml  / 

■i  I Cf  b t-  I OrN 

If  L6  s 0,  the  program  branches  to  8100,  since  no  pivots 
have  been  removed.  Otherwise,  a pivot  has  been  removed  by 
deleting  a row,  so  the  other  elements  of  the  old  pivot 
column  are  tested  for  potential  pivots,  and  if  one  is  found, 
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the  program  branches  to  8110  to  carry  out  the  pivot  operation. 
Otherwise,  the  subroutine  is  finished. 


-i  >10  K 1 » 0 < L C*L  f >=K 
•i  l^f  L )=Lfc 

H 130  53  = 1 
•i  WiC  l u 3uh  * 1 f f 

iitC  »-t  TUM'i 

This  sequence  identifies  the  pivot  row  and  column  as  K1  and 
Ll,  respectively,  records  the  pivot  indices,  sets  the  sensor 
S3  at  +1  to  indicate  a pivot  (as  opposed  to  an  "unpivot, " 
which  would  be  indicated  by  -1) , goes  to  SUB  G160  to  make 
the  pivot,  and  returns  to  the  control  program. 

J . Auxiliary  Subroutines 
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This  subroutine  (through  8390)  controls  the  pivot  operation, 
addressing  SUB  8525  for  the  actual  computation  of  new  values. 
This  sequence  locates  the  first  element  of  the  pivot  row  in 
the  pivot  region  and  in  the  a region,  and  locates  the  pivot 
element.  Then  it  prepares  for  computation  of  the  pivot  row, 
sets  the  multiplier  Q so  that  SUB  8532  will  have  the  effect 
of  dividing  each  element  of  the  pivot  row  by  S3  times  the 
pivot  element,  resets  the  pivot  element  so  that  SUB  8532  will 
have  the  effect  of  replacing  the  original  pivot  element  by 
its  negative  reciprocal,  initializes  15  and  M5  at  the  first 
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elements  of  the  row,  sets  the  stepping  indicators  L7  and 
L8,  and  goes  to  SUB  8532  to  compute. 


BP3C  1 b=.\C 
B £35  Mbs  I < 3t ) 

BP*C  I £=  I 5 + L C I 3 ) * <L  1 - I ) 

B£*b  I-7sLC  \ ?) 
dVAfr  L B =L  ( 13) 

-c  k = c iu  k<  n 

BfbS  U-.SUh  nbPb 
«ytO  I b » J b ♦ j 
BI'VC  JpalB*  1 
3?BC  r*  b = x b ♦ k c i;  > 

BPVt  Nfc.X'1  K 

This  sequence  prepares  for  computation  of  the  rows  in  the 
common  and  6 regions,  initializes  15,  M5 , L7,  and  L8, 
initializes  12  at  the  pivot  column  element  of  the  first  row 
and  loops  on  the  rows,  going  to  SUB  8525  to  compute. 

B3CC  I b=  I C 

H3C  I I V- I b+ l ( 3) * (L I - l ) 

* M.B  L 7=  I 
B 3T.3  = J f 3) 

b 3 r ■<  «b=ivr 

B3ff;  k:/s  k = i i,,  - „ 

1 h k=B)  7 H 
i yr  b 1 1 .'ii.**  bst s 
B33(  I b = 1 b ♦ I ( 1) 

< Vib  .vS:,v;,^  (13) 

B 3B.f  I ?•  = J k + 1 f 1 ) 

3 / 1‘  \i  r.  /,  i ^ 

B J B C M 5 I M ♦ u '( 

*3*0  f-iTUf  -; 

This  sequence  prepares  for  computation  of  the  rows  in  the 
a and  pivot  regions,  except  the  pivot  row,  which  has  already 
been  computed.  It  maintains  the  count  of  the  current  number 
of  explicit  pivots  in  (S)  by  replacing  P8  by  P8  + S3. 
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This  computational  subroutine  is  addressed  from  a number 
of  places  in  the  program.  It  loops  simultaneously  across 
two  parallel  rows  in  (S) , subtracting  from  the  elements  of 
one  row  the  product  of  a multiplier  Q and  the  elements  of 
the  other  row.  In  some  instances,  the  first  two  statements 
of  the  subroutine  are  by-passed.  When  used,  they  initialize 
uhe  multiplier  and  the  pivot-column  element  S(I2)  of  the 
row  being  modified,  so  that  statement  8552  will  have  the 
effect  of  dividing  this  element  by  S3  times  the  original 
pivot  element. 


i 5 7 C U= C 
~i  j i 1 = 15 

-i  59T  H.r  l.=  1 b I.a 

K Jh  l.  (L(*l. i -I  Km  •*».'<( 
<7  t PC.  ■-  = <,*.(  I ) 

I = I ♦Li 
H / A ( \ (•  r J t. 

■1  7 5 C r r.  1 1 it  m 


This  subroutine  is  addressed  for  a specified  row  of  either 
the  <5  region  or  the  pivot  region  of  (S) . It  sums  the  elements 
in  pivot  columns  of  that  row.  It  is  used  only  when  replacing 
a base  row. 
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This  subroutine  reads  one  wing  of  (I)  into  (K)  and  the  other 
into  (L) . When  the  sign  of  SI  is  reversed,  the  wings  are 
interchanged  and  the  problem  is  effectively  transposed. 
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This  subroutine  names  as  simple  variables  a number  of  fre- 
quently used  indexes  associated  with  a target  group  G.  KO, 

LO  are  row,  column  locators  for  (0) . K4,  L4  are  the  number 

of  rows,  columns  in  the  pivot  region  of  (S) . 10  is  the 

location  of  the  first  element  in  the  pivot  region  of  (S) . 

MO  is  the  location  of  the  first  element  in  the  a region  of 
(S) . NO  is  the  location  of  the  first  element  in  the  6 region 
of  (S).  1(1)  is  the  interval  between  elements  of  successive 

rows  in  the  same  column  of  the  pivot  region.  1(3)  is  the 
interval  between  elements  of  successive  columns  in  the  same 
row  of  the  pivot  region. 

^ i 3 u,  .ic=*ck>*k< v>*<7 

r,  1*K  CKC7>*fi>*(  ) 

*r>l.  .J«  = K ( K > + XC  9)  *(  r,  | ♦.•CKf  7 ) ♦ 14.1-v  ) 

-S  'S  *1  ('  i\r  Joi  * 
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This  subroutine  computes  several  indexes  to  (A) . JO  locates 
the  test  control  row  for  group  G.  J1  locates  the  base 
strategy  row  for  group  G.  J9  locates  the  base  strategy  row 
for  the  hypothetical  group  G9. 

K.  Strategy  Subroutine 

This  subroutine  is  regularly  entered  from  statements  1520 
and  1720  in  the  control  program.  There  is  also  a special 
short-cut  entry  from  statement  1100,  used  to  cause  a strategy 
printout  without  changing  resources. 

The  first  portion  of  the  subroutine  sets  up  the  third 
stage  matrix  (U)  and  inverts  it  into  (W> . The  remainder 
carries  the  back  solutions  through  the  third,  second,  and 
first  stages.  The  end  product  is  two  vectors,  (X)  and  (Y) , 
representing,  respectively,  the  first  and  second  columns  of 
the  solution  strategy  expressed  in  terms  of  the  particular 
parameter  being  used  at  the  moment.  (X)  and  (Y)  are  normally 
not  combined  into  a single  solution  strategy  vector  except 
during  printout  in  SUB  5000. 

The  subroutine  serves  eight  cases,  corresponding  to  the 
combinations  of  values  that  SI,  S2,  and  P 2 may  have  in  the 
control  program. 
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This  sequence  zeros  (U) , (W) , and  (Z)  at  the  greater  dimen- 
sion of  (S)  less  the  number  of  pivots  in  (S) . 
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This  sequence  enters  the  common  region  of  (S)  into  (U) . The 
use  of  (I)  indices  means  that  {U)  is  set  up  without  trans- 
position. 
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This  sequence  enters  into  (U)  the  other  elements  of  (S)  at 
the  intersection  of  nonpivot  rows  and  nonpivot  columns.  The 
pivot  index  (0)  controls  the  selection.  Elements  of  the  a 
region  are  entered  by  statement  9150,  of  the  6 reyion  by 
9178,  and  of  the  pivot  region  by  9192. 
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If  S2  = 1,  (U)  is  a nonsingular  matrix  and  may  be  inverted. 

If  S2  = -1,  the  last  row  or  column  of  (U!  consists  of  zeros, 
and  the  deficiency  is  made  up  by  entering  P2  in  the  R5 
element  of  the  last  row  or  column. 
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The  solutions  are  sensitive  to  the  orientation  of  (S)  and 
(W) , so  subroutine  8690  is  called  to  provide  (K)  and  (L) 
indices,  and  (Z)  is  set  as  either  (W)  or  the  transpose  of 

(W)  . 
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This  sequence  computes  the  third-stage  solution  from  (Z)  and 
the  right-side  elements  of  (S) , storing  the  solution  in  (C) 
and  (F). 

-.-a!.;  k »•  - = l i.  *.  ( ) ) 

i 'K  ( *.  t * ) - i.  ( « > 

-i  *i  1 1 r(«  )jp(h) 

< ",  i t i - 

This  sequence  enters  the  marginal  returns  from  (C)  and  (F) 
directly  into  (X)  and  (V) . 

1 1 -•„(•>♦  i c * ) • • 
l i -•(!>♦  1 

i ) 

r . : ~ i I . i 1 

This  sequence  begins  the  back  solution  by  a loop  on  G 
extending  to  line  9900.  For  G = 1,  it  initializes  Jl  to 
locate  the  base  strategy  in  (A)  and  II  to  locate  the  inter- 
cept elements  in  (X)  and  (v) . It  then  initializes  K2  to 
locate  the  last  elements  of  (C)  and  (F)  already  transferred 
to  (X)  and  (Y).  For  each  group  in  turn,  it  goes  to  SUB  8760 
to  get  specific  indexes. 
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This  sequence  begins  the  back  solution  for  elements  of  (X) 
and  (Y)  corresponding  to  excess  strategies.  The  loop  on  I 
extends  to  line  9700.  If  there  is  a pivot  in  row  Kl  of  (S) , 
the  program  branches  to  9610;  otherwise,  it  continues. 


■/-I.  f(  1 > s r ( - ; > 


If  the  Kl  row  of  (s)  is  not  a pivot  row,  the  corresponding 
solution  elements  (X)  and  { Y } are  entered  directly  from  (C) 
and  (F)  . 
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If  the  Kl  row  of  (S)  is  a pivot  row,  the  corresponding 
elements  of  (X)  and  (Y)  must  be  computed  from  (C) , (F) , 

and  all  those  elements  of  (S>  in  the  column  of  that  pivot 
which  are  not  themselves  in  any  pivot  rov:.  The  computation 
runs  first  over  the  marginal-value  elements  of  (C)  and  (F) . 
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The  computation  is  completed  by  running  over  other  elements 
of  (C)  and  (F)  pertaining  to  target  group  G. 
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The  back  solution  for  intercept  value,  X(I1)  and  Y(I1),  and 

for  base  strategy,  X(J)  and  Y(J),  begins  with  their  initialization. 
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This  sequence  completes  the  back  solution  for  intercept  and 
base  strategy,  using  previously  computed  elements  of  (X)  and 
(Y)  and  stored  elements  of  (R)  and  (A) . 
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This  sequence  steps  indexes  for  the  next  target  group.  After 
completing  (X)  and  (Y) , it  returns  to  the  control  program. 

The  END  statement  is  not  a part  of  the  operating  program 
but  is  required  by  the  computer  system. 


V.  PATH87A  COMPUTER  PROGRAM 


One-sided  optimization  problems  can  be  run  on  the  PATH87 
program  by  specifying  no  defended  target  groups  and  no  defense 
weapon  types,  i.e.,  G4 , Dl  = 0.  However,  much  of  the  capa- 
bility of  the  program  is  not  needed  for  these  simpler  problems 
and  can  be  cut  out  to  save  running  time. 

The  PATH87A  program  is  designed  for  one-sided  optimizations. 
Being  derived  from  PATH 8 7 , it  has  a similar  general  structure 
but  is  less  than  half  the  length.  Excess  capabilities  have 
been  cut  out  and  many  procedures  simplified. 

A . Simplifying  Ideas 

By  definition,  the  defense  solution  strategy  is  a null 
stL^tegy,  which  need  not  be  computed  or  stored.  Since  there 
are  no  excess  elementary  defense  strategies,  the  pivot 
regions  cf  (R)  and  (S)  are  nonexistent,  and  the  second-stage 
solution  procedures  can  be  eliminated,  along  with  the  pivot 
index  (0) . If  the  components  of  the  attack  solution  strategy 
are  redefined  as  numbers  of  targets,  rather  than  as  fractions 
of  a group  of  targets,  the  scaling  factor  T(G)  can  be 
eliminated  from  (S) , which  reduces  round-off  error  as  well 
as  computational  time.  The  one-sided  nature  of  the  problem 
results  in  most  of  the  subroutines  having  to  work  only  one 
way  rather  than  two.  Consequently,  the  reversible  indexing 
system  is  not  needed,  and  we  can  eliminate  (K)  and  (L) . 

B . Program  Listing 

The  following  listing  has  explanatory  comments  added 
where  appropriate: 
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5 0 I/AJA  b*  3 
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7 C L-AJA  * 8#  . 5,  . b,  » 7,  . t , . 3#  . /#  . t , . a,  .9,  . s,  . 3,  . 4,  .p,  . S 

Only  two  dimensions,  G1  and  Al,  need  to  be  specified.  Line  80 
is  omitted. 


icr  Mr  Mirn.K  ico.w  icr>.M<prp>,  upc?>,<  < ic  i>,.K?r?-> 

ICb  01*  f(  15), l-C  lb), ef  1 b >, i. < |S,  |S), '/<!*>,  |b) 
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This  sequence  omits  matrices  (D) , (K) , (L) , (0)  , (R) , and 

(Z)  . 


(G)  is  reduced  in  size  since  only  the  a region  of  (S) 
needs  indexes. 

(Q)  is  enlarged  to  store  the  value  associated  with  every 
elementary  attack  strategy,  not  just  the  base  strategies. 

(M)  is  used  to  locate  the  base  value  for  group  G in  (Q)  . 

(N)  is  reduced  to  store  only  the  number  of  rows  for 
each  target  group  in  (S) . 

(J)  is  added  as  a locator  for  the  attack  region  of  (A) . 
It  is  divided  into  two  regions,  J(l)  to  J(G9)  locating  base- 
strategy  rows  and  J(G9+1)  to  J(2*G9)  locating  control  rows 
for  each  group.  The  defense  region  of  (A)  is  not  used. 
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This  sequence  reads  dimensions  of  the  problem  and  initializes 
indices.  In  actuality,  the  program  could  be  improved  by 
replacing  the  (I)  elements  by  simple  variables  throughout. 


»(  ( 

l :■ 

IM 

* • 

i 3 

V 

j 1C 

i ■ .. 

f i 

= 

1 

1 

1 1 

VI 

rl 

I 

< 

1 ) 

» 

V < ( 

) 

.ViC 

I h 

I \ i 

1 

< ( 

) 
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This  sequence  reads  target  data  and  initializes  the  (G)  , 
(M) , and  (J)  indices. 
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1 

( 

> i 1 N 1 
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y.  7C 

X 1 

c 
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VIM 

C 

\M.  1 

( 

This  sequence  reads  probabilities  of  kill  and  stores 
probabilities  of  survival  in  (P) . 
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1CCC  FM.Ml  ••F-l^^a?" 

ICIC  INPUT  Pl#I-3 

to  pc  if  f i<»  c thfn  ic*o 

1030  blUP 

Meaningful  inputs  for  PI  are  positive  integers  to  denote 
attack  weapon  types,  0 to  terminate  the  run,  and  -1  to  g^L 
a printout  of  marginal  values.  Their  use  is  illustrated  in 
Section  VI. A. 


1C4C  d|sSUN<M) 

1 0 b f.  IF  h 1 = - 1 1 H LfM  1 1 C 0 
1 C7C  I * < PI)  ) 

I C7  t,  v V.  1 » < 1 - *-£•  > 

i (."i c u=F«io:i-'(i'i)) 

KVC  If  F*»  0 1 HK.m  1130 

nor  uiiuF  4pst 

II  I ( Gi;  1 G 1 4 ( C 
1130  K!K  Fo=1  7/  4P 

ll^C  IF  nc  I <P3>  + f-b>  = H 7mf\  r/oo 
1130  IF  (■(  \ ( r 3 ) * I-  b ) > P I THIN  I I H C 
1 1 1 r \M7  b 5 
1170  1 

1 1*0  F P P=  f-'c  IP  hb  ilH-1 

1 f-b4  I = 1 ♦ F- 

\?bt  hub  K = c 70  PP 

1PF  0 > C 1*1)  = .->(  1 ) 
i ?f.p  m i ) = o 
1 ^ e 3 i = i + ic:h) 

1 76  4 .MLX7  K 
1 ? 7 0 \/Fx7  ^ 

lPHC  F OF  . I = 1 <•  0-  3 > 7 k.  J<  f,9  ) - I ( 76  ) b 7 F F \ ( 7 6 h 

I 79  C FGI-  L-.I+  I < 3 1 ) 1;:  J*hb  b7F.  P -1 

1 300  ^(L+Ds^d) 

13  10  N M i L 
I 3F0  (•(.)*  hb  > = 0* 

1330  M F X T ,1 

1 340  PC  1 < ?3>*Ib)  = Pl 

1 350  £?=£?♦ 1 

134  0 mnsi(.7i)+i 

This  sequence  initializes  a path  segment  in  the  same  general 
fashion  as  in  PATH87,  but  with  many  simplifications.  The 
number  of  attack  weapon  types  actually  in  play  is  denoted  by 


A2  rather  than  by  1(4) 


1 5(  i. 

. 1 - - 1 

1 t.;-  i 
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l * .» C 
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* :.‘- 

: r " 
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: l *■  i 

i -if  ( 

j .,;*.!*• 

/(  ( r 

J -i  1 f 

-li. 

l * 

The  1500 

and  1700  s 

in/'*. 

« < i i ) 

s ^ i 

Kit 

r f l ( ) * .•  s 

1 -f  <■  1. 

•i  = i ^ 

* « i ; 

1 -5  - < 

i *• 

\ ( * r 

> - ; < 


This  sequence  terminates  the  path  segment. 


i-  C V.  I r | .*  1 i-.r.  . 

/■(->(  U:( 

r > f I = I I , I I 
> l V . , I f 1. 1 r • I * < ) 

tif'..  . I 1 - . I ( i ) 
y : n < i = - . r f < i t ) 
t < * f i i - 1' 

; (;•< : ^ I - i 7.:  * •- 

r If  I < I = ’ i - ' ( I > « (-  ( . I ) + | ) 

; i ( j f i ~t  i - f ( i >*/-(-. 'i*i  > 

.III  ^A|  | 

This  sequence  begins  the  test  control  section  for  a defense 
strategy  by  initializing  XI  and  Yl.  The  procedure  permits 
SUB  9000  to  be  greatly  simplified. 
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f t P r ( K JU"  ?p  (■( 

P I3C  l.ClC  P19C 
p'  i ^ c Uwi,uh  a^rc 
p i sr  <-i  ik  pr/r 
put  *=i 

ir  I V f IK  /•(  > )*  A(  JC»w>  Hh  PP3C 

PHC  f.-bui-  vrr 
P i*t  ic  r i>  -.  c f i i -c t ^ :-;>c 

P c C.  L uisu1"  4 C f. C 

PPIC  Uk.-L'J^-  33C  f 

ppr-r  ik  #•{  i o i > j<  = *-<  jc*  i c 3 1 > > imkn  fun 
ff  ■»(  IK  /<<-:  HH  PUC 

f f 4*J  | <4  = I ^ ♦ I-  l 

p : ‘.j  c m k ^ i ( . 

f f t ( ^ t t U /•  Ni 

This  sequence  completes  the  test  control  section  for  a 
defense  strategy  with  essentially  the  same  flow  illustrated 
in  Figure  11.  Three  exceptions  should  be  noted.  SUB  3000 
is  omitted  since  the  defense  strategy  is  unchangeably  a null 
strategy.  SUE  4000  is  addressed  directly  instead  of  through 
SUB  3700.  SUB  3820  is  replaced  by  a single  statement  stepping 
14. 

f i><  r h'h  l-rj  If  (-J 
PblC  rubU^  3CCC 
P t 7 C MX]  0 
p t ",  f Kill  /•  .'M 

The  test  control  section  for  an  attack  strategy  calls  on 
SUB  3000  only. 
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The  3000 

-series  subroutines 

these  are  simplified. 
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The  simplified  value  function  computes  X2  directly.  The 
variable  Y2  is  omitted  since  it  would  always  be  zero. 
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bGGC  Ofc>  fna< I ) *M*y ( l ) ♦ a c I > 

STIC  V*G 

bflb  I V bid  THtN  ‘iWC 

bCPC  huh  V-*l  Tl  Af- 

bCtC  hhlNl  A(l(?3)*v.>i-X(k) 

bC7C  >N  F X T w 
SCVb  KtluKN 

The  print  subroutine  begins  with  one  branch  printing  only 
the  part  of  a defense  r trategy  that  specifies  marginal  values 
for  the  attack  weapons. 
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The  other  branch  prints  an  attack  strategy  with  no  marginal 
values  or  intercept  values.  VS  is  computed  as  the  sum  of 
the  products  (numbers  of  targets  times  expected  surviving 
value  per  target).  Illustrative  printouts  are  in  Section  VI. A. 
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The  basis  subroutines,  6000  and  7000,  realize  the  greatest 
saving  of  statement  lines  by  the  elimination  of  the  second 
solution  stage  and  all  the  complicated  algorithms  associated 
with  it. 
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The  matrix  (U)  is  set  up  from  only  one  region  of  (S)  instead 
of  four.  If  Si  = -1,  the  program  goes  to  9411  to  compute  a 
defense  strategy.  Otherwise,  it  continues  at  9330. 
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This  sequence  begins  to  compute 

an  attack  strategy  by  storing 

the  third 

-stage  solutions  in  <C) 

and  (F).  The  program 

continues 

at  9490. 

j 

*•*-»! 

f i.  t-  \ 1 

J 

*•<.13 

s 

u U 

' ( ( 

**  1 *> 

vh  «•  i •< 

» M H 

* = J ( 3 1 ) 

* 4?  C 

► . r L = 1 i l (■' 

4 «•  «o 

■*  s y ♦ I ( 3 k ) 

/ 4.  4.  f 

hi  r m = ) h < r 

A(4):^(K)».S(''),,.(4,U 

1 

y i.4«4 

\ t X 1 K 

1 

•# 

N t x I L 

1 

* *4-  S 

- ^ lur\ 

1 

This  sequence  computes  a defense 

strategy,  marginal 

values  1 

only,  storing  the  third-stage  solutions  directly  in 

(X)  ard 

(Y) ; then 

it  returns  to  the  control  program. 
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This  sequence  finishes  computing  an  attack  strategy,  less 
marginal  and  intercept  values;  then  it  returns  to  the  control 
program. 
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VI . EXAMPLES 


In  this  section  we  will  apply  the  PATH87  and  PATH87A 
programs  to  solve  the  problem  of  allocating  weapons  to  point 
targets  in  a variety  of  cases  selected  to  illustrate  features 
of  the  solution  method  and  typical  properties  of  solutions. 

A.  One-Sided  Optimizations;  PATH87A 

The  first  four  cases  are  one-sided  optimizations  involving 
the  allocation  of  attack  weapons  to  undefended  targets. 

Program  PATH87A  is  preferred  for  this  type  of  problem, 
although  PATH 8 7 may  be  used. 

The  dimensions  of  each  case  are  given  by  the  number  of 
target  types  (Gl)  and  the  number  of  attack-weapon  types  (Al) . 
The  cases  are  arranged  in  order  of  increasing  complexity. 

Case  1;  (Gl,Al)  = (1,1) 

The  technique  can  be  illustrated  by  running  a simple 
case  with  interaction  between  computer  and  operator.  To  begin, 
we  assume  that  the  source  program  is  stored  on  disk. 

L i AU  1-A1HB7A 
h E ADY 

The  operator  commands  that  the  source  program  be  loaded  into 
core,  and  the  computer  responds  when  it  has  done  this  and 
is  ready  for  further  commands . 

bG  CAT  A 1,1 
60  OAT A 1 CC, l 
VC  DATA  .3 

r.  UM 

The  operator  modifies  the  source  program  and  commands  that 
it  be  run.  In  this  example,  the  data  statements  specify  1 
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target  type,  1 attack-weapon  type,  100  targets,  unit  value 
per  target,  and  0.8  single  shot  kill  probability. 

PATH87  A 

li  1*  A I = I 1 

T,  Vs 
100  1 
tKT~ 

.8 

The  computer  prints  a heading  and  begins  the  run  by  printing 
the  data  specifications  for  the  record. 

M,P3?  f 

? i,5c  '( 


After  initializing,  the  computer  calls  for  input  of  PI  and 
P3  to  specify  a path  segment.  By  typing  "1,50",  the 
operator  directs  that  attack  weapon  type  1 be  changed  from 
its  present  quantity  (initially  0}  to  50. 


1 

C 50  50 

1 5C  10. 

60. 


At  the  end  oi  the  path  segment,  the  computer  prints  a 
solution,  four  lines  in  this  case.  The  first  line  identifies 
the  target  type.  The  second  line  says  that  the  elementary 
strategy  a = 0 (i.e.,  no  weapons  arc  assigned)  is  used  on 
each  of  50  targets  and  that  the  expected  surviving  value  of 
these  is  50.  The  third  line  says  that  the  elementary 
strategy  a = 1 (i.e.,  one  weapon  is  assigned  to  each  target) 
is  used  on  each  of  50  targets  and  that  the  expected  sur- 
viving value  of  these  in  10.  The  last  line  gives  the  total 
expected  surviving  value. 


“ * H 


MjK'i.' 

? - 1*  r 

» 

The  computer  calls  for  a new  path  segment,  and  the  operator 
types  "-1,0”,  a special  input  used  to  direct  a printout  of 
marginal  values  in  PATH87A.  The  computer  says  that  the 
marginal  value  for  attack  weapon  tyne  1 is  -.8. 
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In  the  same  way  as  before,  the  number  of  weapons  is  increased 
to  150,  with  a new  solution  strategy  and  a new  marginal  value 
of  -.16. 
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The  force  size  is  increased  to  250,  and  the  marginal  value 

becomes  -.032. 
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A "0,0"  input  terminates  the  run,  with  the  computer  printing 
the  total  central  processing  time  used. 
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The  results  are  illustrated  graphically  by  the  drawdown 
curve  of  Figure  13.  Here  the  resource  space  is  the  A-axis, 
divided  into  regions  of  length  100  by  regional  boundaries  at 
the  points  0,  100,  200,  etc.  Elementary  solution  strategies 
throughout  the  first  region  are  a - 0 and  a - 1,  throughout 
the  second  a - 1 and  a = 2,  etc.  The  solution  surface  is 
the  drawdown  curve  itself,  a straight  line  over  each  region. 
The  tangent  line,  with  its  intercept  and  slope  (or  marginal 
value) , is  uniquely  defined  over  the  interior  of  each  region 
but  is  indeterminate  at  regional  boundaries,  where  it  may 
have  any  slope  between  those  of  the  two  adjacent  regions. 


Solution  Surface.  VS 


FIGURE  13 

SOLUTION  OF  CASE  1 


Case  2;  (G1,A1)  * (1,2) 

A case  with  two  attack  weapon  types  illustrates  some 
additional  features.  For  the  following  runs,  we  assume  that 
PATH87A  has  already  been  loaded  into  core. 


b(’  DA TA  1,P 
( C DATA  1 DC/  1 
1 C DATA  • i , . t 


The  operator  commands  a run  with  1 target  type,  2 attack 
weapon  types,  100  targets,  unit  value  per  target,  and  0.8 
and  0.6  single  shot  kill  probabilities. 
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The  computer  prints  the  heading  and  data  record. 
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For  the  first  path  segment,  the  input  ”1,150"  produces  a 
solution  identical  with  one  of  the  Case  1 solutions. 


r U A3? 
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The  second  path  segment  adds  weapon  type  2.  In  the  solution, 
there  are  now  three  elementary  strategies,  given  by  the 
vectors  (2,0),  (1,2),  and  (2,1). 

bl#f3? 

? -I,C 

1 CCCC-(f--Cr 

V 4CUU-CP 

The  marginal  values  for  both  attack  weapons  are  obtained  by 
the  special  input  *'-l,0". 

» l,  K3? 

v t 

Cf  u-or.Cli:  C.  b 

The  run  ends  with  a time  printout. 

Another  run  with  the  same  data  statements  can  be  made 
simply  by  commanding  RUN. 


/-  PIH8  /£ 

(n  1 , £ 1 = I r 

1 » v- 

ICC  I 

hKl  = 

.3  . t 

> i » r3V 
? JPC 

) 

l np. 

P 'ri-  'Ur 

3b. P 

In  this  example,  the  order  of  weapon  input  is  reversed, 
attack  weapon  type  2 being  raised  to  120  on  the  first  path 
segment. 
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On  the  second  path  segment,  attack  weapon  type  1 is  added. 

In  the  printout  of  elementary  strategies,  the  weapon  types 

are  arranged  in  their  original  numerical  order  regardless  of 
the  path  sequence.  The  solution  strategy  at  the  point  (150, 
120)  is  identical  with  that  of  the  preceding  run. 

In  this  case,  the  resource  space  is  the  (At ,A2 ) -plane. 

It  is  divided  into  triangular  regions  as  shown  by  the  solid 
lines  on  the  regional  map  of  Figure  14.  The  solution  surface 
must  be  imagined  in  a third  dimension,  but  the  values  are 
shown  in  parentheses  at  points  of  integral  density,  correspond- 
ing to  corners  of  the  regions.  The  solution  surface  over 
each  region  is  a plane.  The  elementary  solution  strategy 
at  the  corner  of  a region  is  a uniform  strategy,  the  same  on 
all  targets.  On  a boundary,  the  solution  is  a mix  of  two 
corner  strategies.  In  a regional  interior,  the  solution  is 
a mix  of  three  corner  strategies. 

The  pattern  of  regions  is  associated  with  the  kill 
probabilities  of  the  problem.  In  this  case,  the  pattern  was 
determined  by  making  a run  along  the  dashed  exploratory  path 
in  the  figure.  The  exploratory  path  is  one  that  goes  back 
and  forth,  like  ploughing  a field.  Ke  want  a printout  at 
every  critical  point  of  the  path,  that  is  at  every  regional 
boundary  and  at  the  end  of  each  path  segment.  For  this 
purpose,  the  program  must  be  modified. 
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H7C 
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The  operator  adds  lines  1745  and  1750  to  get  a printout  at 
every  critical  point,  and  deletes  1870  to  delete  the  usual 
printout  at  the  end  of  a path  segment. 
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The  solution  at  (50,0) , the  end  of  the  first  path  segment, 
is  identical  with  that  in  Case  1. 
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On  the  second  path  segment,  running  from  (50,0)  to  (50,280), 
the  computer  prints  a solution  strategy  at  the  first  regional 
boundary  encountered.  The  three  elementary  strategies  locate 
the  corner  points  of  a region.  Since  one  of  the  strategies, 
i.e.  , (0,0),  applies  to  zero  targets,  the  path  is  crossing 
the  opposite  regional  boundary,  specifically  at  the  point 
(50,50) , where  VS  = 30. 
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The  (0,0)  elementary  strategy  has  been  deleted  and  (0,2) 
is  added.  The  path  has  reached  the  point  (50,100),  where 
it  is  again  on  a regional  boundary,  and  the  (0,1)  elementary 
strategy  is  about  to  be  deleted. 


V ~=  1.  'I'lV't  i 

The  path  continues,  encountering  regional  boundaries  at 
(50,150),  (50,200),  and  (50,225).  Here  a small  roundoff 
error  appears  in  the  printout  of  VS,  which  should  be  read 
as  5.8 
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At  the  boundary  (50,250),  a number  of  small  roundoff  errors 
appear.  The  correct  reading  is  obviously: 
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The  program  continues  along 

the 

dashed  line 
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The  path  segment  ends  at  an  interior  point  of  a region.  No 
roundoff  errors  appear  in  the  printout. 
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The  third  path  segment  traverses  four  regions  before 
reaching  a boundary  at  (113j,280>.  Neither  the  force  level 
nor  the  strategy  is  physically  feasible,  although  they  are 
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mathematically  correct  in  the  path  method  of  solution 
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The  path  segment  continues  to  its  end  at  the  interior  point 
(150,280).  Here  the  strategy  is  physically  feasible  when 
roundoff  errors  in  the  printout  are  overlooked. 
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On  the  fourth 

path 

segment. 

the 

so  the  value  surviving  is  increasing.  At  (150,120),  the 
end  of  the  segment,  the  solution  is  identical,  except  for 
arrangement,  with  those  of  the  earlier  runs  of  this  case. 
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The  tilth  path  segment  continues  in  the  direction  of  the 
fourth,  with  VS  still  increasing. 
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The  force  increases  again  on  the  sixth  path  segment. 
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The  force  continues  to  increase  on  the  seventh  path  segment, 

and  VS  reaches  its  lowest  level  of  the  entire  run  at  the  end 
point  (250,120). 
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The  last  path  segment  reaches  point  (150,120)  from  the  'ight 
in  the  map  of  Figure  14,  with  the  same  solution  as  from  each 
of  the  other  directions. 


f 1 > K3  t 

■'  c,i 

<’  h u - 3 b Co  : t ♦ (• 

The  exploratory  run  is  terminated  with  sufficient  data  to 
construct  the  map  of  Figure  14.  By  inspecting  the  figure, 
we  can  see  that  the  path  has  traversed  each  region  at  least 
once  and  some  regions  as  many  as  four  times. 

Although  the  path  method  theoretically  finds  the  same 
solution  regardless  of  the  direction  of  approach  to  a point, 
the  program  PATH87A  does  not  always  do  so,  because  of 
arbitrary  limitations  of  the  "floating  lid"  test-control 
process  discussed  in  Section  TV.  The  discrepancy  is  illustrated 
by  the  two  runs  that  follow. 
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In  the  first  run,  the  correct  VS  = 5.44  at  (150,60)  is 
obtained  by  inputting  (2,60)  and  then  (1,150). 
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In  the  second  run,  with  the  input  order  reversed,  the  answer 
VS  = 5.52  is  not  optimal,  because  the  elementary  strategy  (0,3) 
was  barred  from  testing  when  the  boundary  at  (150,50)  was 
crossed,  since  the  "floating  lid"  would  not  permit  testing 
either  weapon  at  a level  higher  than  1 above  its  previous 
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When  the  same  run  is  continued  by  withdrawing,  input  (2,20), 
and  then  readvancing,  input  (2,60),  the  correct  solution  is 
obtained,  because  the  lid  has  floated  high  enough  to  permit 
the  elementary  strategy  (0,3)  to  be  tested  the  second  time 
the  boundary  is  crossed.  The  "floating  lid"  process  needs 
improvement  in  order  to  achieve  absolute  precision,  as  well 
as  to  save  time. 

Case  3 • (G1,A1)  ■ (1,5) 

Some  additional  features  can  be  illustrated  by  a case 
with  five  attack-weapon  types. 
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The  resource  space  is  now  three  dimensional  and  is  divided 
into  tetrahedral  regions  with  four  corners  and  triangular 
faces.  The  solution  contains  four  elementary  strategies. 
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Each  new  weapon  type  adds  a dimension  to  the  resource  space 


and  an  elementary  strategy  to  the  solution.  The  general  rule 
for  one-sided  optimizations  is  that  the  number  of  elementary 


strategies 

is 
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to  the 
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of  object  types 

plus  the 
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The  final  resource  space  has  five  dimensions  and  is  divided 
into  six-cornered  regions  called  simplexes.  Roundoff  errors 
in  the  printout  of  the  target  column  can  be  corrected  to  58, 
12,  7,  17,  2,  and  4,  giving  an  integral  solution  strategy. 
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The  marginal  values,  or  slopes  of  the  tangent  hyperplane, 
are  all  very  small  at  this  attack  level. 

We  cannot  guarantee  that  the  final  solution  strategy  is 
absolutely  the  best  integral  strategy,  since  the  "floating 
lid"  test  process  may  have  arbitrarily  excluded  a better 
elementary  strategy  than  one  of  those  selected.  However,  we 
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can  find  a lower  bound  for  VS  by  supposing  that  the  aggregate 
kill  potential  of  the  force  could  be  spread  uniformly  over 
all  the  targets.  This  lower  bound  is  VS  = 5.4589,  computed 
by  hand.  Clearly  the  integral  strategy  of  the  computer  run 
couldn't  be  improved  very  much,  if  at  all. 

Case  4:  (G1,A1)  = (5,3) 

Multidimensionality  may  extend  to  both  target  and  weapon 
types,  as  in  the  following  run. 
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The  five  target  groups  consist  of  1 target  of  value  10,  15 
targets  of  value  8,  etc.,  a total  of  96  targets  of  aggregate 
value  390.  The  3 weapon  types  have  different  kill  probabilities 
against  targets  of  the  different  groups.  The  first  weapon 
type  is  generally  the  best,  but  the  third  weapon  type  is  best 
against  group  5 targets. 
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On  the  first  path  segment,  50  type  3 weapons  are  used.  There 
are  two  elementary  strategies  on  target  group  4,  but  only  one 
elementary  strategy  on  each  of  the  other  groups.  In  the  sense 
of  Figure  13,  that  means  there  is  an  interior  strategy  on 
group  4 and  corner  strategies  on  the  other  groups. 
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The  marginal  value  is  determined  by  the  interior  strategy 
on  group  4. 
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On  the  second  path  segment,  50  type  2 weapons  are  added  to 
the  force,  and  some  reallocation  of  type  3 weapons  takes 
place.  In  the  sense  of  Figure  14,  there  are  corner  strategies 
on  groups  1,  3,  and  4,  edge  strategies  on  groups  2 and  5,  and 
no  interior  strategies. 
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The  marginal  value  for  weapon  type  2 is  determined  by  the 
strategy  on  target  group  2,  and  that  for  weapon  type  3 by 
the  strategy  on  target  group  5. 
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Now  50  type  1 weapons  are  added,  and  the  other  weapons  are 
reallocated.  We  note  that  all  the  type  3 weapons  are  used 
against  target  group  5,  where  they  are  more  effective  than 
types  1 and  2.  There  are  corner  strategies  on  two  groups 
and  edge  strategies  on  three.  Boundary  strategies  of  some 
sort  are  typical  in  multigroup  cases. 

The  strategies  on  groups  3 and  4 are  not  physically 
feasible,  since  they  call  for  the  use  of  an  elementary 
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strategy  on  a fractional  number  of  targets.  If  desired,  the 
solution  can  be  adjusted  by  hand  to  a feasible  one,  for 
instance: 


2 0 0 

0 3 0 
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2 0 0 

VS-  67.672, 


3 1.62 

6 2.304 

1 .288 
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which  is  close  to  the  printed  value.  Ordinarily,  the 
discrepancy  can  be  ignored,  since  it  is  a local  one,  as  is 
shown  by  the  following  path  segment. 


r 1> 

! a i 
1 


V 

r 

i 

. & 

i 

t J 

r 

(• 

t>.  1 t 

'c 

( 

7 

) . (-'r 

C 

,7 

7 . 

r\  (.  ■. 

1 

>• 

t . 

1 'jj  • 

f . S 

t 1 

r 

v 

7 . 

. 7 b 

c 

i 

i 

tc 

bC 

Vo=  <7.  .**15 


If  one  type  2 weapon  is  added  to  the  force,  the  strategy 
becomes  physically  feasible.  The  changes  taking  place  show 
that  physically  feasible  strategies  will  result  from  48,  51, 
54,  etc.,  type  2 weapons  and  infeasible  ones  from  49,  50,  52, 
53,  etc.  The  user  will  have  to  decide  whether  extra  precision 
is  worth  the  trouble. 


124 


h 1 1 K V! 

v -»#r 

I -.Ab 

3 * 

h ij  f-3? 

? f 

IlMj-jM!:  i J-.f' 

The  marginal  values  show  *hat,  at  the  final  force  level 
(50,51,50),  one  extra  type  3 weapon  would  add  more  to  force 
effectiveness  than  one  of  either  of  the  other  types. 

B.  Two-Sided  Games;  PATHS 7 

The  next  four  cases  are  two-sided  games  involving  the 
allocation  of  defense  weapons  as  well  as  of  attack  weapons. 
Program  PATH 8 7 must  be  used  for  these. 

In  addition  to  the  dimensions  Gl  and  Al,  we  must  also 
specify  the  number  of  defended-target  types  (G4)  and  the 
number  of  defense-weapon  types  (Dl) . 

Case  5:  (G1,G4,A1,D1)  = (1,1, 1,1);  Perfect  Weapons 


This  perfect-weapon  case  is  analyzed  in  Reference  1, 
where  explicit  solution  formulas  are  given.  It  is  included 
here  to  illustrate  the  numerical  solution  by  the  PATH87 
program. 

Si  i./ i A ] * i » l » i 

(.<  C V i-  i b 1 C < * 1 
vr  i >r  ir-  i 

HC  "fit  \ 

:■  l:-: 

The  first  line  of  data  input  specifies  the  four  dimensions. 
The  added  fourth  line  specifies  the  probability  of  intercept 
by  the  defense  weapon  type.  The  other  lines  are  the  same  for 
both  PATHS 7 and  PATH87A. 
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After  printing  a heading  and  recording  data,  the  computer 
calls  for  the  operator  to  input  a path  segment.  Some  attack 
must  be  specified  before  any  defense,  so  the  operator  inputs 
50  type  1 attack  weapons. 
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The  form  of  the  strategy  printout  differs  in  two  respects 
from  that  of  PATH87A.  Here  the  first  line  shows  an  intercept 
value  of  50,  and  the  elementary  strategy  lines  omit  the 
expected-survival  column  of  PATH87A. 
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The  operator  inputs  175  type  1 defense  weapons  (the  - sign 
indicates  defense) . 
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The  computer  prints  the  defense  solution  strategy  at  the 
end  of  the  second  path  segment)  that  is,  at  the  point 
(50jl75).  The  first  line  says  attack  weapon  type  1 has  a 
marginal  value  of  -.225 — that  is,  the  slope  of  the  tangent 
plane  in  the  direction  of  the  attack  coordinate.  The  second 
line  says  that  on  target  group  1 the  intercept  value  is 
100-— that  is,  the  ordinate  [i.e.,  value  at  the  point 
(0)175)]  of  the  plane  tangent  to  the  solution  surface  at 
(50)175).  The  next  five  lines  give  elementary  defense 
strategies  and  the  number  of  targets  for  each.  The  last 
line  gives  the  expected  value  surviving  at  the  point  (50)175) 
as  88.75,  which  is  computed  as  the  intercept  value  plus  the 
marginal  value  multiplied  by  the  number  of  attack  weapons. 

The  weapon  allocation  itself  exhibits  the  equal  steps, 
i.e.,  22.5  targets,  characteristic  of  perfect  weapon  cases. 
When  these  fractional  allocations  are  viewed  as  mixed 
strategies,  they  are  physically  feasible,  since  numerous 
convex  combinations  of  integral  allocations  can  be  found 
that  are  precisely  equivalent.  Perhaps  the  simplest  combi- 
nation consists  of  the  strategy 
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with  a frequency  of  .5  and  the  strategy 
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with  a frequency  of  .5. 
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The  operator  calls  for  an  attach  strategy  printout  by 
inputting  an  attack  weapon  at  no  change  from  its  current 
force  level. 
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Although  the  attack  force  level  is  unchanged,  the  attack 
solution  strategy  has  changed.  Before,  there  were  50  targets 
attacked.  Now  tiwt  a defense  has  been  introduced,  the  attack 
is  concentrated  in  an  equal  step  s.^ategy  on  only  20  targets. 
Notice  that  the  marginal  value  per  defense  weapon  is  .05, 
the  intercept  value  at  (50;0)  is  80,  and  VS  = 80  + 175  x .05  = 
88.75,  as  it  should. 

A complete  solution  at  the  point  (50; 175)  is  provided  by 
this  attack  strategy  and  the  preceding  defense  strategy.  The 
point  happens  to  lie  in  a region  that  we  call  defense  dominant, 
where  some  targets  are  not  attacked. 
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A complete  solution  at  the  point  <275; 175)  is  obtained  by 
first  increasing  the  attack  force  to  275  and  then  calling 
for  a defense  printout.  This  point  lies  in  a region  of 
attack  dominance,  since  every  target  is  attacked  by  at 
least  one  weapon.  The  roundoff  arrors  are  annoying,  but 
they  can  be  easily  adjusted  in  this  case,  for  instance  by 
reading  both  22.5006  and  22.4999  as  22.5. 
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The  run  is  terminated  in  the  usual  fashion. 

The  same  problem  can  be  run  with  slightly  more  generality 
as  a density  problem  by  normalizing  inputs  as  in  the  follow- 
ing run. 
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The  data  change  specifies  1 target  instead  of  100. 
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PATHS  7 
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The  first  path  input  is  a density  of  attack  of  .5  instead 
of  50  weapons.  Subsequent  inputs  will  be  scaled  the  same 
way. 
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The  results  of  the  density  run  should  be  compared  item  by 
item  with  those  of  the  preceding  run.  The  solution  strategies 
and  values  can  be  scaled  to  apply  to  any  number  of  targets 
of  any  unit  value.  An  extra  benefit  of  the  density  technique 
is  the  reduction  of  roundoff  errors. 


Case  6: 


(G1,G4,A1,D1)  = (1, 1/1,1);  Imperfect  Weapons 
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The  solution  strategies  at  the  point  (17/3; 11/3)  are  identical 
with  the  canonical  strategies  pictured  in  Reference  2 
(Figure  23,  p.  172)  except  for  notational  differences. 

PATH87  can  be  used  for  exploratory  runs,  as  in  the  follow- 
ing example,  which  also  illustrates  certain  noteworthy 
features  of  the  solution  process. 
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The  program  must  ba  modified  to  print  a complete  spectrum 
of  strategies  along  the  path. 
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The  first  path  segment  runs  from  the  origin  to  the  point 

(.1/0) . 
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The  first  strategy  is  a defense  strategy  at  the  origin , in 
the  usual  form  except  for  the  third  line,  which  is  an 
intentionally  incomplete  representation  of  the  null  alloca- 
tion at  the  origin.  The  only  important  line  is  the  one  that 
specifies  the  initial  marginal  value  for  the  attack  as  -.8. 

. lrc 
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. 1 

Next  is  an  attack  strategy,  the  rule  being  an  alternation  of 
attack  and  defense  strategies  in  the  spectrum  along  a path 
segment.  As  no  defense  weapon  has  yet  been  introduced,  the 
printout  of  marginal  value  is  suppressed  and  the  first  print 
line  shows  an  intercept  value  of  .92,  the  same  as  VS.  Since 
the  strategy  allocates  all  prescribed  weapons,  the  path 
segment  is  terminated. 
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The  input  for  the  next  path  segment  calls  for  raising  the 
defense  force  to  a density  of  1 weapon  per  target.  Whenever 
a weapon  type  is  introduced,  the  first  strategy  is  an 
opposition  strategy,  so  here  an  attack  strategy  is  printed. 
It  is  the  same  as  the  strategy  terminating  the  preceding 
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segment  except  that  it  now  contains  a marginal  value  of  .06 
for  the  weapon  type  just  introduced. 
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Next  is  a defense  strategy  in  standard  form.  Since  it  occurs 
at  a force  level  of  .592593,  which  is  before  the  end  of  the 

path  segment,  we  know  that  a regional  boundary  has  been 

encountered  at  this  force  level.  As  both  elementary  strate- 
gies, 6=0  and  6=1,  are  shown  to  apply  to  more  than  zero 
targets,  neither  can  be  deleted;  so  we  infer  that  a new 
elementary  attack  strategy  must  be  added  to  the  Q-basis  at 
the  boundary.  The  value  of  the  game  at  the  point  (.1;. 592593) 
is  .955556. 
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The  printout  of  the  attack  strategy  on  the  boundary  shows 
that  the  elementary  strategy,  a = 2 , has  been  added,  and 
the  attack  is  now  concentrated  on  a smaller  fraction  of  the 
target  system.  The  VS  printout  here  should  be  disregarded, 
because  the  program  computes  VS  using  the  defense  force  level 
recorded  in  the  ($)  matrix,  and  this  matrix  is  not  updated 
until  the  end  of  the  path  segment.  Actually,  this  attack 
strategy  is  a solution  strategy  for  a range  of  defense  force 
levels,  beginning  at  .592593,  where  VS  = .944615  + (.592593)* 
(.0184614)  = .955555,  the  same  as  printed  for  the  preceding 
defense  strategy.  The  range  includes  the  following  defense 
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force  level*  which  we  can  anticipate  by  computing  VS  = .944615 
+ (1.0)  (.0184614)  * .963076. 
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The  path  segment  ends  with  a defense  strategy  that  includes 
the  elementary  strategy  6=2.  A complete  solution  at  the 
end  point  (.1;1)  consists  of  this  defense  strategy,  the 
preceding  attack  strategy,  and  the  value  of  the  game,  which 
is  computed  here  as  VS  = 1 + (. 1 )(-. 369231)  = .963077. 

The  course  of  events  along  the  path  segment  can  be 
rationalized  in  the  following  terms: 

First— at  the  beginning,  all  of  the  attack  weapons 
are  on  separate  targets,  since  there  is  no  defense. 

Second — As  defense  weapons  are  added,  they  are 
allocated  to  separate  targets  with  an  expected  saving  of 
.06  targets  per  weapon. 

Third — When  the  defended  fraction  of  targets 
reaches  .592593,  it  becomes  just  as  good  for  the  attack  to 
put  2 weapons  as  1 on  a target. 

Fourth — At  this  point,  the  attack  shifts  to  a 
strategy  that  kills  just  as  many  targets  but  offers  the 
defense  the  smallest  possible  saving  for  any  additional 
weapons,  i.e.,  .0184614.  The  new  strategy  is  so  balanced 
between  the  targets  attacked  with  1 weapon  and  those 
attacked  with  2 weapons  that  the  defense  can  do  equally 
well  by  defending  1 target  with  2 weapons  as  by  defending 
2 targets  with  1 weapon  each. 
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Fifth — The  defense  strategy  does  not  change  at  this 
point,  but  as  new  defense  weapons  are  added  to  the  force, 
the  balance  between  the  elementary  strategies  6=1  and 
6 ■ 2 is  maintained  so  that  the  attack  is  unable  to  profit 
by  changing  allocation. 

It  should  be  recognized  that  each  strategy  printed  out 
along  the  path  segment  is  a solution  strategy  when  paired 
with  the  preceding  strategy  of  the  opposition,  the  following 
strategy  of  the  opposition,  or  any  convex  combination  of  the 
preceding  and  following  strategies.  Each  printed  attack 
strategy  is  valid  for  some  range  of  this  path  segment,  and 
convex  combinations  of  them  represent  transitions  at  a 
boundary  point.  Each  printed  defense  strategy  is  valid 
only  at  a point  of  the  path  segment,  and  convex  combinations 
of  them  represent  solutions  at  intermediate  points.  In 
this  sense,  the  printed  strategies  constitute  a complete 
spectrum  of  solutions  along  the  path  segment. 
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The  path  input  calls  for  the  attack  density  to  be  increased 
from  .1  to  2.6.  Since  no  new  weapon  type  is  introduced,  the 
first  strategy  printed  is  an  attack  strategy  at  the  critical 
point  (1.625;1),  the  first  boundary  encountered  on  the  path 
segment.  Here  all  targets  come  under  attack,  and  the 
elementary  strategy  a = 0 leaves  the  Q-basis.  The  path  is 
leaving  a region  of  defense  dominance  an  i entering  a region 
of  attack  dominance. 
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At  the  boundary,  the  defense  makes  a transition  to  a strategy 
for  the  next  region.  The  new  strategy  concentrates  on 
defending  fewer  targets  and  offers  a numerically-reduced 
marginal  value  to  the  attack.  Since  no  elementary  strategy 
is  zeroed  out,  we  infer  that  the  next  attack  strategy  will 
contain  a new  elementary  strategy. 
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The  next  attack  strategy  includes  a = 3,  as  inferred.  The 
attack  force  is  now  at  a critical  value  of  2.25.  No  elementary 
strategy  is  zeroed  out,  so  we  infer  that  the  next  defense  will 
have  an  added  elementary  strategy. 
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The  defense  does  have  four  elementary  strategies,  but  6-1 
is  immediately  zeroed  out  during  the  transition  at  the 
boundary.  The  defense  is  now  more  concentrated,  and  the 
marginal  value  is  reduced  in  magnitude. 
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The  path  reaches  another  boundary  at  the  critical  attack 
size  of  2. 480768 , where  a * 1 is  zeroed  out  and  every  target 
is  attacked  by  at  least  two  weapons. 
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The  defense  transitions  to  a more  concentrated  strategy, 
again  reducing  the  magnitude  of  the  marginal  value. 
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At  (2.6;1),  the  end  point  of  the  path  segment,  a complete 
solution  is  given  by  this  attack  strategy  and  the  preceding 
defense  strategy.  The  VS  can  be  computed  from  either  one, 

VS  = .502606  + (2.6)  (-.11826) 

= .0229565  + (1)  (.172173) 

= .19513  . 

This  is  the  same  value  given  in  Reference  2 (page  159)  . The 
spectrum  of  strategies  on  the  preceding  path  segment  is 
identical,  except  for  notation,  with  part  of  the  spectrum 
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pictured  in  Reference  2 (Figure  18,  page  15?) . For  the 
next  path  segment,  we  will  now  run  a portion  of  the  other 
spectrum  appearing  in  Reference  2 (Figure  19,  page  158). 
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The  path  input  calls  for  increasing  the  defense  to  a density 
of  3.5  weapons  per  target.  The  first  defense  strategy  is 
at  the  critical  density  of  1.4. 
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A new  attack  and  a new  defense  have  four  elementary  strategies 
each. 
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Here  the  attack  transitions  to  a more  concentrated  strategy, 
leaving  some  targets  unattacked.  The  path  is  entering  a 
defense  dominant  region. 
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This  critical  defense  strategy  displays  a new  phenomenon. 
When  the  defense  force  has  increased  to  a size  of  3.0365, 
it  becomes  worthwhile  to  defend  every  target  with  at  least 
one  weapon,  so  the  6=0  line  leaves  the  Q-basis. 
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The  attack  responds  to  the  new  situation  by  attacking  every 
target  at  least  once,  and  reducing  the  marginal  value  for 
new  defense  weapons,  which  can  only  be  added  above  the  1 
level. 
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The  defense  increases  to  a point  where  it  pays  the  attack 
to  use  5 weapons  on  some  targets. 
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The  attack  transitions  to  a strategy  where  more  targets  are 
attacked  by  1 weapon  only,  the  balance  of  the  force  being 
concentrated  on  fewer  targets. 
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The  path  segment  ends  at  (2. 6; 3. 5),  an  interior  point  of  a 
region.  The  run  is  terminated  after  2.1  seconds  of  central 
processing  time,  including  1.3  seconds  of  time  to  compile 
the  program  before  running. 

Figure  15  shows  a regional  map  for  these  weapon  types. 

The  general  pattern  of  interlocking  rectangular  regions  is 
characteristic  for  one  imperfect  attack-weapon  type  and  one 
imperfect  defense-weapon  type  although  the  precise  locations 
of  the  boundaries  vary  with  the  probabilities  of  kill  and 
intercept.  The  pattern  is  the  same  as  that  shown  in  Figure  17 
of  Reference  2.  The  solution  surface  over  each  rectangular 
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The  defense  increases  to  a point  where  it  pays  the  attack 
to  use  5 weapons  on  some  targets. 
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The  attack  transitions  to  a strategy  where  more  targets  are 
attacked  by  1 weapon  only,  the  balance  of  the  force  being 
concentrated  on  fewer  targets. 
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The  path  segment  ends  at  (2.6;3.5),  an  interior  poinc  of  a 
region.  The  run  is  terminated  after  2.1  seconds  of  central 
processing  time,  including  1.3  seconds  of  time  to  compile 
the  program  before  running. 

Figure  15  shows  a regional  map  for  these  weapon  types. 

The  general  pattern  of  interlocking  rectangular  regions  is 
characteristic  for  one  imperfect  attack-weapon  type  and  one 
imperfect  defense-weapon  type  although  the  precise  locations 
of  the  boundaries  vary  with  the  probabilities  of  kill  and 
intercept.  The  pattern  is  the  same  as  that  shown  in  Figure  17 
of  Reference  2.  The  solution  surface  over  each  rectangular 
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region  is  a hyperbolic  paraboloid  with  rulings  in  the 
coordinate  directions,  as  discussed  in  Reference  2.  These 
features  contrast  with  the  triangular  regions  and  plaaar 
solution  surfaces  already  shown  for  cases  with  two  attack- 
weapon  types. 

Case  7:  {Gl,G4,Al,Dl)  » (1,1,2, 2) 

The  following  density  run  illustrates  the  multiweapon 
case. 
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The  attack  weapons  are  conventionally  arranged  in  order  of 
increasing  PKT,  0.6  and  0.8;  the  defense  weapons  in  order 
of  decreasing  PI,  0.9  and  0.75. 
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The  first  two  path  segments  input  attack  weapon  type  2 and 
defense  weapon  type  2.  As  these  are  the  same  types  used  in 
Case  6,  the  solution  at  the  point  (G,.L;0,1)  is  also  the 
same  as  in  the  exploratory  run  of  that  case. 
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At  the  point  (0,.1;1,1)  the  defense  solution  strategy  is 
stepped  in  form,  with  some  targets  being  undefended  and 
some  defended  by  from  1 to  4 weapons.  The  bettor  inter- 
ceptor is  concentrated  on  about  half  the  targets  and  the 
poorer  interceptor  spread  out  over  about  86  percent  of 
the  targets. 
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At  the  point  (.6,. 6; 1,1),  the  attack  strategy  shows  the 
poorer  weapon  concentrated  on  about  one  fourth  of  the  targets 
and  the  better  one  spread  out  over  about  44  percent  of  the 
targets.  The  levels  of  attack  range  from  0 to  4 weapons  per 
target. 
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A printout  of  the  defense  strategy  shows  only  minor  changes 
resulting  from  the  increased  attack. 
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The  run  is  terminated. 
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The  resource  space  is  four  dimensional,  so  we  can't  show 
a regional  map  of  it.  Actually,  we  don't  know  enough  to 
discuss  the  regions  in  much  detail.  The  space  is  bounded 
by  six  coordinate  planes:  a map  of  the  (At ,A2) -plane  has 
triangular  regions  like  Figure  14  with  the  axes  reversed;  a 
map  of  the  (A2 ,D2 ) -plane  is  the  same  as  Figure  15;  maps  of 
the  (Aj ,Dj )-,  (A»,Da)-,  and  (A2 ,Dj ) -planes  have  rectangular 
regions  similar  to  those  of  Figure  15;  and  a map  of  the 
(Dj ,D2) -plane  cannot  be  defined  since  VS  = 1,  regardless  of 
the  defense  strategy,  when  there  is  no  attack.  Intuitively, 
we  expect  the  interior  regions  to  have  both  rectangular  and 
triangular  faces,  and  this  is  confirmed  by  a few  exploratory 
paths,  which  have  also  shown  that  some  of  the  multidimensional 
regions  overlap. 

Case  8:  (G1 ,G4 ,A1 ,Dl)  = (10,4,1,1) 

. ■ ■ ■ { 
The  multiple-target-type  case  is  illustrated  by  the 

following  run. 
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There  are  10  target  groups,  only  four  of  which  are  defended 
(by  convention,  the  first  four).  Each  group  has  10  targets, 
the  total  number  of  targets  is  100,  and  the  aggregate  value 
is  550.  The  single  attack  weapon  type  has  a different 
probability  of  kill  for  each  group  of  targets. 
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On  the  first  path  segment,  the  attack  lays  down  165  weapons 
as  if  all  the  targets  were  undefended.  The  marginal  value 
to  the  attack  has  not  been  printed,  but  it  can  easily  be 
computed  from  the  strategy  on  target  group  2 as  (9)(l-.6)1 
(-.6)  = - .864. 
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On  the  second  path  segment,  the  defense  allocates  55  weapons. 
VS  increases  from  59.0426  to  144.487.  The  marginal  value  to 
the  attack  increases  to  -1.25112,  as  determined  by  the  defense 
strategy  on  the  first  four  groups. 
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A printout  of  the  attack  strategy  at  the  point  (165; 55) 
shows  target  group  10  no  longer  attacked  and  target  group  5 
under  reduced  attack.  The  shift  of  20  weapons  from  these 
two  groups  to  the  defended  groups  takes  advantage  of  the 
higher  marginal  values  provided  by  the  defense.  Some  shift 
of  the  attack  from  undefended  to  defended  targets  is  usual 
as  the  defense  begins  to  build  up.  If  the  defense  gets  very 
strong,  however,  the  marginal  values  may  decrease,  and  some 
or  all  of  the  attack  weapons  may  shift  back  to  the  undefended 


targets.  The  operator  of  the  program  should  be  alert  for 
a shift  of  all  the  attack  weapons  to  undefended  targets , 
since  the  resulting  indeterminate  defense  could  lead  to  a 
failure  of  the  program.  The  program  ought  to  have  a 
procedure  for  terminating  the  path  segment  in  such  an  event, 
but  the  current  version  does  not. 
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The  run  terminates. 
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VII.  OTHER  APPLICATIONS 


The  illustrative  cases  in  Section  VI  involved  no 
substantial  changes  to  the  PATH 8 7 and  PATH87A  programs, 

since  only  data  input  and,  in  a few  examples,  printout  were 
affected. 

Other  types  of  resource  allocation  problems  can  be  run 
if  the  programs  are  modified  to  fit  the  problem.  Some  that 
have  actually  been  run  include  such  features  as: 

-1  Terminal  defenses 

~2  Mix  of  area  and  terminal  defenses 

-3  Overlapping  area  defenses 

-4  Sensitivity  to  errors  in  force  estimation 

-5  Decoy  weapons 

-6  Decoy  targets 

-7  Weapon  allocations  in  sets  of  2 or  more 

-8  Area  targets 

-9  Area-mobile  targets 

“10  Targets  of  variable  value 

“11  Time-phased  attacks 

-12  Defense-suppression  attacks 

-13  Transportation  cost  matrices 

14  Simultaneous  variation  of  several  resource  types 

-15  A "floating  box"  to  control  testing. 

In  this  section,  we  will  discuss  general  methods  of  adapting 
the  programs  to  various  types  of  problems. 

Preceding  page  blank 


A.  Choice  of  Program 


Generally/  if  the  problem  is  a one-sided  optimization  it 
can  best  be  handled  by  modifying  FATH87A.  Examples  are 
terminal-defense  problems  and  transportation  problems.  Two- 
sided  games  require  PATH87. 

B.  Choice  of  Variables 

Resource  types,  object  types and  measures  of  value  must 
be  chosen  in  an  appropriate  fashion,  keeping  the  dimensions 
of  the  problem  as  small  as  possible.  In  many  cases,  the 
choice  is  obvious;  in  others,  some  judgment  must  be  exercised. 

In  the  overlapping  area-defense  case,  each  set  of  inter- 
ceptors with  the  same  coverage  may  be  considered  a resource 
type,  so  that  there  wj 11  be  as  many  resource  types  as  there 
are  interceptor  sets.  Each  group  of  targets  covered  by  the 
same  interceptor  sets  may  then  be  considered  as  an  object 
type  (but  only  if  these  targets  are  otherwise  identical). 

An  incidence  matrix  may  be  used  to  specify  the  match-up  of 
interceptor  types  and  object  types. 

In  another  type  of  game , the  object  types  were  groups 
of  bomber  bases,  the  defense  resources  were  bombers  to  be 
allocated  to  bases  and  the  attack  resource  types  were  succes- 
sive waves  of  a time -phased  attack  on  the  bases,  with  the 
potential  value  per  base  decreasing  as  surviving  bombers 
were  launched  during  the  attack. 

In  transportation  problems,  the  destinations  are  con- 
sidered as  object  types  and  the  requirement  at  each  destina- 
tion as  the  number  of  objects  of  that  type.  The  sources  are 
considered  as  resource  types  and  the  availability  at  each 
source  as  a quantity  of  resources.  If  preferred,  sources 
and  destinations  can  be  reversed  in  meaning,  with  the  idea 
of  reducing  the  number  of  resource  types. 


C.  Value  Function 

A value  function  for  point  targets  is  built  into  PATHS 7 
and  PATH87A.  It  must  be  replaced  when  it  is  inappropriate. 

For  area  targets,  the  square  root  law  has  been  used  as 
a value  function. 

For  mobile  targets,  the  region  of  mobility  may  be  con- 
sidered as  a single  target  with  initial  value  equal  to  the 
number  of  mobile  targets  in  the  region.  Depending  on 
whether  the  region  is  linear,  areal,  or  spatial,  a value 
function  can  then  be  developed  to  compute  expected-value 
surviving  for  all  allowable  elementary  attack  strategies. 

In  the  bomber-basing  case,  the  initial  value  of  a base 
is  zero,  and  the  maximum  value  surviving  is  the  number  of 
bombers  in  the  elementary  defense  strategy  for  that  base. 

The  value  function  takes  into  account  the  time-phasing  of 
the  elementary  attack  strategy  and  the  bomber  launches. 

In  the  defense-suppression  case,  the  value  function 
incorporates  a suboptimization  of  the  split  between  attack 
weapons  allocated  to  radars  and  those  allocated  to  targets. 

In  the  transportation  problem,  the  value  function  is  a 
precomputed  matrix  of  costs,  one  element  for  each  source- 
destination  pair.  Costs  may  be  in  miles,  dollars,  or  what- 
ever measure  is  to  be  minimized. 

D.  Test  Controls 

"Floating  lid"  test  controls  are  used  in  PATH87  and 
PATH87A.  They  should  be  changed  whenever  better  ones  can 
be  devised  for  a problem. 

In  the  transportation  problem,  there  is  never  more  than 
one  resource  unit  (weapon)  allocated  to  a destination  require 
ment  unit  (target) . Since  tests  of  two  or  more  would  be 
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completely  meaningless,  the  test  controls  can  be  greatly 
simplified. 

In  the  mobile  target  cases  actually  run,  there  were  only 
five  allowable  elementary  strategies  on  targets  of  each 
group.  Again,  a simple  test  control  procedure  was  used. 

Zn  the  bomber-basing  case,  the  results  of  early  runs 
showed  certain  systematic  patterns  in  the  solution  strategies; 
that  knowledge  was  used  to  change  the  test  controls  to  save 
time  on  subsequent  runs . 

E.  Economy  in  Modifications 

If  a special  type  of  problem  is  going  to  be  run  only  a 
few  times,  it  is  usually  economical  to  make  minimal  modifica- 
tions to  one  of  the  basic  programs. 

On  the  other  hand,  if  many  runs  are  planned,  it  may  be 
economical  to  make  more  extensive  modifications  that  take 
advantage  of  simplifying  features  of  the  problem.  For 
example,  if  each  elementary  attack  strategy  is  limited  to  a 
single  weapon  type,  each  one  can  be  described  by  identifying 
the  weapon  type  and  the  quantity,  i.e.,  by  using  two  numbers 
instead  of  a complete  vector.  In  this  case,  the  (A)  matrix 
illustrated  in  Figure  7,  Section  IV,  can  have  its  blocks 
reduced  to  two  elements  instead  of  A1  + 1(27).  Going  even 
further,  if  the  quantity  is  always  one,  as  in  the  transpor- 
tation problem,  then  only  the  resource  identifying  number 
need  be  stored  in  (A).  Of  course,  modifications  of  this 
type  require  that  references  to  (A)  be  modified  throughout 
the  program. 

F . Potential  Improveme nts  and  Applications 

Many  applications  of  the  PATH  programs  have  been  demon- 
strated. As  time  goes  <:n,  we  expect  more  to  be  found.  We 
also  expect  that  substantial  programing  improvements  will  add 
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to  the  efficiency i flexibility,  reliability,  and  usefulness 
of  the  method. 
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The  PATS  Betted  offer*  unique  advantages  of  speed  end 
flexibility  in  solving  problems  facing  defense  analysts,  and 
it  is  hoped  that  publication  of  this  report  through  the 
National  Technical  Information  Service  of  the  Defense 
Documentation  Center  will  make  this  method  sore  widely 
available.  Also , the  method  has  features  which  can  be 
applied  to  auiny  problems  of  resource  allocation  facing 
nondefense  planners. 
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